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We construct two exactly soluble lattice spin models that demonstrate the importance of three- 
loop braiding statistics for the classification of 3 D gapped quantum phases. The two models are 
superficially similar: both are gapped and both support particle-like and loop-like excitations similar 
to that of charges and vortex lines in a Z2 XZ2 gauge theory. Furthermore, in both models the particle 
excitations are bosons, and in both models the particle and loop excitations have the same mutual 
braiding statistics. The difference between the two models is only apparent when one considers the 
recently proposed three-loop braiding process in which one loop is braided around another while 
both are linked to a third loop. We find that the statistical phase associated with this process is 
different in the two models, thus proving that they belong to two distinct phases. An important 
feature of this work is that we derive our results using a concrete approach: we construct string and 
membrane operators that create and move the particle and loop excitations and then we extract the 
braiding statistics from the commutation algebra of these operators. 


I. INTRODUCTION 

The discovery of quantum Hall states and more re¬ 
cently, topological insulators^, has taught us that there 
are many different types of gapped quantum many body 
systems. In order to understand the relationship between 
these systems, it is useful to divide them into classes in 
such a way that the members of each class share the same 
qualitative properties. Typically these classes are defined 
as follows: two gapped Hamiltonians H and H' are as¬ 
signed to the same class if they can be adiabatically con¬ 
nected to one another — that is, if there exists a one pa¬ 
rameter family of interpolating Hamiltonians H(s) with 
(1) H (0) = H , H( 1) = H' and (2) a finite energy gap for 
all s, 0 < s < 1. The different classes of Hamiltonians 
are then called “gapped phases.” 

The precise definition of gapped phases depends on 
what type of systems we wish to consider. For example, if 
we are interested in systems with particular symmetries, 
then it is natural to assign H and H' to the same phase if 
there exists an interpolating Hamiltonian H(s ) which is 
both gapped and invariant under the relevant symmetries. 
Including such symmetry constraints typically leads to a 
finer classification of gapped phases, as illustrated by the 
example of topological and conventional insulators^ 

In this paper, we consider the coarsest possible classi¬ 
fication scheme: that is, we do not impose any symmetry 
constraints and we say that two gapped Hamiltonians 
H and H' belong to the same phase if they can be adi¬ 
abatically connected by any interpolating Hamiltonian 
H(s) with local interactions. Our starting point is a ba¬ 
sic question: how can we tell whether or not two gapped 
Hamiltonians belong to the same phase? 

This question has an appealing answer in the case of 
two dimensional (2D) systems^ To determine whether 
two gapped 2D Hamiltonians belong to the same phase, 
one can simply compare the braiding statistics^ of their 
quasiparticle excitations. If the braiding statistics data 
do not match, then the two Hamiltonians must belong to 
different phases since braiding statistics cannot change 


under an adiabatic deformation. Conversely, if the braid¬ 
ing statistics do match, then we can almost conclude that 
the two Hamiltonians belong to the same phase. To reach 
this conclusion, we simply need to compare one other 
quantity, namely the thermal Hall conductance^. Indeed, 
according to a plausible (but unproven) conjecture, if two 
Hamiltonians have the same braiding statistics and the 
same thermal Hall conductance, then they must be adia¬ 
batically connected to one another in the absence of any 
symmetry constraints. 

In the 3D case, our understanding is much more lim¬ 
ited. One way to attack the classification problem is 
to try to generalize the concept of quasiparticle braid¬ 
ing statistics to the 3D case. The simplest generalization 
begins with the observation that many 3D Hamiltoni¬ 
ans support loop-like excitations in addition to particles. 
Given this observation, we can consider several differ¬ 
ent types of braiding statistics. First, we can look at 
the exchange statistics of particle-like excitations. These 
exchange statistics can take only one of two values for 
each particle: every particle must be either bosonic or 
fermionic^ Second, we can consider the statistical phase 
associated with braiding a particle around a loop^""— 
Finally, we can look at the statistical phase associated 
with braiding one loop excitation around another (Fig. 
[Tt)JA--ia if we combine all of these types of braiding 
statistics, we can indeed distinguish many different 3D 
gapped phases. 

Interestingly, however, this data is incomplete: Refs. 
M and argued that we also need to consider the sta¬ 
tistical phase associated with a three-loop braiding pro¬ 
cess. In this process, one loop is braided around another 
loop while both are linked to a third loop (Fig. (Ud). It is 
unclear whether the three-loop braiding data is the last 
piece of the puzzle or whether there exist further distinc¬ 
tions between 3D gapped phases that can only be seen 
if we consider other braiding processes or other types of 
probes. However, either way, three-loop braiding statis¬ 
tics has already proven to be useful in certain cases 

One weakness of previous studies of three-loop braiding 
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statistics is that this quantity has only been calculated 
using indirect and abstract arguments. For example, 
Ref. E3 computed the three-loop braiding statistics of 
3D Dijkgraaf-Witten gauge theories using a dimensional 
reduction argument which relates the braiding statistics 
of loops in 3D Dijkgraaf-Witten models to the braiding 
statistics of particles in 2D Dijkgraaf-Witten models^. 
The approaches of Refs. and k> were also indirect: 
Refs. [HI and Qjl computed loop braiding statistics by 
relating this quantity to modular transformations on a 
three dimensional torus. 

In this paper, we study three-loop braiding statistics 
using a more concrete approach. We focus on two exactly 
soluble lattice models and we compute their three-loop 
braiding statistics by explicitly implementing the loop 
braiding process on the lattice. Our approach is simi¬ 
lar to how quasiparticle braiding statistics is commonly 
calculated in 2D lattice models^— : we construct mem¬ 
brane operators that create and move loop-like excita¬ 
tions, and then we extract the three-loop statistics from 
the commutation algebra of these operators. 

The two spin models that we analyze provide an 
explicit demonstration of the importance of three-loop 
braiding statistics for distinguishing 3D gapped phases. 
Indeed, we show that the models share the same parti¬ 
cle exchange statistics as well as the same particle-loop 
and loop-loop braiding statistics. The only difference be¬ 
tween the models is that they have different three-loop 
braiding statistics. Thus, it is only this quantity that 
reveals that the two models belong to different phases. 

The models that we study are not completely new and 
have appeared previously in the literature in different 
forms. In particular, the first model is essentially iden¬ 
tical to the 3D generalized toric code mode l 21 i 23 ~ — cor¬ 
responding to the group Z 2 x Z 2 . Thus, the low energy 
properties of this model are similar to that of conven¬ 
tional Z 2 x Z 2 gauge theory2£. The second model can 
be thought of as a different type of Z 2 x Z 2 gauge the¬ 
ory. More specifically, this model can be obtained by 
starting with the spin model in Ref. H3 which describes 
a nontrivial symmetry-protected topological phased - — 
with Z 2 x Z 2 symmetry, and then coupling this system 
to a Z 2 x Z 2 lattice gauge field2&. We will explain this 
connection in more detail in a separate publication^. In 
addition to its connection with the spin model in Ref. [27l, 
we believe that the second model belongs to the same 
phase as one of the exactly soluble Z 2 x Z 2 Dijkgraaf- 
Witten models i 2Q i 36 This conjecture is based on the fact 
that the braiding statistics in the two systems seem to 
match one another. 

The connection between our models and symmetry- 
protected topological phases is not accidental: we specif¬ 
ically designed our models to be equivalent to gauged 
symmetry-protected topological phases, because accord¬ 
ing to the results of Ref. Q3 and [RJ, such gauge theories 
can support different types of three-loop braiding statis¬ 
tics. The reason that we chose to gauge the spin model 
from Ref. |27] is that this is one of the simplest known 


(a) (b) 
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FIG. 1. (a) Two-loop braiding process, (b) Three-loop braid¬ 
ing process. The gray curves show the paths of two points on 
the moving loop. 



FIG. 2. Both models are built out of two species of spins. 
The blue spins a v live on the plaquettes p of the cubic lattice 
while the red spins ap live on the plaquettes p of the dual 
cubic lattice. 

models for a 3D symmetry-protected topological phase 
with unitary symmetry group. 

The rest of the paper is organized as follows. In section 
mi we introduce the two exactly soluble 3D spin models 
that we will analyze. In section ITTT1 we study the particle¬ 
like and loop-like excitations of these models by explic¬ 
itly constructing the string and membrane operators that 
create and move these excitations. We then compute the 
braiding statistics of these particles and loops in section 
EH Some technical details can be found in the appen¬ 
dices. 


II. MODELS 

A. Hilbert space for the models 

The models that we will discuss are spin-1/2 systems 
made up of two species of spins: “blue spins” and “red 
spins.” The blue spins live on the plaquettes p of the 
cubic lattice and will be denoted by <j p while the red 
spins live on the plaquettes p of the dual cubic lattice 
and will be denoted by ap (Fig. [2]). In this notation, the 
S z eigenstates |{cr|, crj}) provide a complete basis for the 
Hilbert space. 

We will often find it convenient to describe spin states 
using an alternative language based on “membranes.” In 
the membrane language, each S z eigenstate |{<Tp,crJ}) 
corresponds to a spatial configuration of red and blue 
membranes on the cubic lattice and dual cubic lattice. 
The dictionary between spin states and membrane con¬ 
figurations is as follows: if a z = —1 we say that the 
plaquette p is occupied by a blue membrane, while if 
a z = +1 we say that the plaquette p is empty. Similarly, 
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if crj = — 1 then p is occupied by a red membrane while if 
cr| = +1, then p is empty. In this way, each spin state can 
be equivalently described as a membrane state. We will 
label our membrane states as \Xb,X r ) where X 5 denotes 
the subset of plaquettes that are occupied by blue mem¬ 
branes, and X r denotes the subset of plaquettes occupied 
by red membranes. 

B. Ground state wave functions 

The two models that we will discuss have been engi¬ 
neered to have particular ground states. These ground 
states are easiest to describe if we assume an infinite 
(non-periodic) geometry. In such a geometry, the ground 
state of the first model is 

|*o)= £ \x b , Xr) (1) 

closed Xb,X r 

where the sum runs over all closed membrane states 
\Xb,X r ). Here, by a closed membrane state, we mean 
a membrane configuration (Xb, X r ) in which all the blue 
and red membranes form closed surfaces, i.e. surfaces 
without boundaries. More precisely, a closed membrane 
state is defined to be a state in which every edge l in the 
cubic lattice and every edge l in the dual cubic lattice is 
adjacent to an even number of occupied plaquettes. 

The ground state of the second model has a similar 
form 

|*i>= £ (-l) N ^ x »’ x ^\X b ,X r ) (2) 

closed Xb,X r 

where again the sum runs over all possible closed mem¬ 
brane states \Xb,X r ). The quantity N g (Xb, X r ) is de¬ 
fined as follows: for each closed membrane configura¬ 
tion (Xb,X r ), the intersections between the red and 
blue membranes X 5 D X r form a collection of discon¬ 
nected closed curves, which we will call “green loops.” 
The quantity N g (Xb,X r ) is defined to be the number of 
the green loops in X 5 D X r . As discussed in the intro¬ 
duction, the ground state |Ti) is closely related to the 
ground state of the spin model in Ref. [23 which describes 
a nontrivial symmetry-protected topological phase with 
Z 2 x Z 2 symmetry. 

C. Sheared cubic lattice 

The reader may notice that there is a technical prob¬ 
lem with the above definition of N g {Xb , X r ): the problem 
is that the closed membrane condition allows for mem¬ 
brane configurations in which an edge l is adjacent to four 
occupied plaquettes. Geometrically such configurations 
correspond to the case where two blue membranes touch 
one another along the edge l. This membrane touching 
is problematic because it means that the “green loops” 
defined by the intersections of red and blue membranes 



FIG. 3. (a) The sheared square lattice is formed by shifting 
the corners of the square plaquettes (the blue dots) according 
to (i,j) —> [i + ej, j) with e > 0. (b) Top view of the sheared 
cubic lattice. The solid/dashed squares denote cubes on two 
neighboring layers. The sheared cubic lattice is formed by 
shifting the corners of the cubes according to (z, j, k) —>• {i + 
ej + ek,j + e'k, k) with e' > e' > 0. 


can also touch one another at corners. As a result, there 
is some ambiguity in determining the number of discon¬ 
nected green loops corresponding to a membrane config¬ 
uration (Xb,X r ). 

To deal with this issue, we now describe a way to 
infinitesimally deform the cubic lattice so as to elimi¬ 
nate membrane touching. Before describing this defor¬ 
mation, we first warm up with an analogous deforma¬ 
tion of the square lattice. The basic idea is to think of 
the square plaquettes that make up the square lattice as 
rigid blocks that can be shifted around. The deforma¬ 
tion we have in mind corresponds to shifting the position 
of the square plaquettes so that their corners move from 
(i,j) (i + ej, j). The resulting lattice, which we call 
the sheared square lattice is shown in Fig. [3^a). 

We are now ready to consider the cubic lattice. In 
this case, we think of the cubes that make up the cubic 
lattice as rigid blocks and we shift these cubes so that 
their corners move from 

(i, j, k) (i - 1 - ej + e'k, j + e'k, k ) (3) 

where e' > e > 0. The resulting sheared cubic lattice is 
shown in Fig. [3^b). 

Now let us imagine performing the same shearing de¬ 
formation to both the original cubic lattice and the dual 
cubic lattice. We can then deform an arbitrary closed 
membrane configuration (Xb,X r ) on the cubic lattice 
(and dual cubic lattice) to a membrane configuration on 
the sheared cubic lattice (and sheared dual cubic lattice). 
The result is a closed membrane configuration without 
any membrane touching of any kind. In what follows, we 
will always think of membrane configurations as living 
on the sheared cubic lattice rather than the cubic lattice 
so we don’t have to worry about membrane touching. 
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D. Hamiltonians 

1. Ho 

We are now ready to write down the Hamiltonians for 
the two models. The Hamiltonian for the first membrane 
model is a sum of four terms, 

= - £> - £4- - £ B c - (4) 

1 i c c 

where the indices /, l run over the links of the cubic lattice 
and dual cubic lattice respectively, while c, c run over the 
“cubes” of the cubic lattice and dual cubic lattice. To 
define the operators Ai, B c , B&, it suffices to explain 
how they act on the membrane basis states |X&, X r ). The 
Ai , Af operators are given by 

Ai = -(i + s&i) , Aj- = -(1 + (5) 

where 

M\x b ,x r ) = (~l) N '\X b ,X r ), 
£/ t \X b ,Xr)-**(-l) N '\X b ,X r ). ( 6 ) 

Here Ni and Nj are defined to be the number of occupied 
plaquettes adjacent to l and l respectively. The B®,B® 
operators are defined by 

B c ° = i(l + B c °) , Bl= l -{l + Bl) (7) 

where 

B° c \X b ,X r ) = \X b + c,X r ), 

Bl\X b ,X r ) = \X b ,X r + c). (8) 


(a) (b) (c) (d) 



FIG. 4. (a) A blue cube intersects a red membrane; their 

intersection consists of a single red loop, (b) A blue cube 
intersects/overlaps with a blue membrane; their intersection 
consists of a single blue region, (c) A blue cube intersects with 
both a red membrane and a blue membrane; their intersection 
consists of one red loop and one blue region, (d) For the 
membrane configuration shown in (c), the integer m c — 1 
because there is one red loop while n c = 2 because there are 
two intersections (denoted by green dots) between the red 
loop and the boundary of the blue region (thick blue line). 


where these products run over the six plaquettes adjacent 
to c, c respectively. 


2. H 1 

The Hamiltonian for the second membrane model is 
similar: 

^ = (H) 

l l c c 

The first two terms are defined as above, while the last 
two terms are defined by 


Here the notation X^ + c is meant to denote a kind of Z 2 
addition on membrane configurations. More specifically, 
given the membrane configuration X 5 , the configuration 
X\) + c is obtained by flipping the occupation numbers of 
the six plaquettes of the cube c: that is, one changes the 
plaquettes from unoccupied to occupied and vice versa. 

There is a simple physical picture for the A and B 
in terms of the membrane language: the A terms fa¬ 
vor closed membrane configurations, while the B terms 
provide an amplitude for the membranes to fluctuate. 
Together these terms ensure that the ground state is a 
superposition of many different closed membrane config¬ 
urations. 

Alternatively, we can express the A and B operators 
in the spin language: 

Ai = \d + U^) ’ A i=\( 1 + 1>J) (9) 

pel 

where these products run over the four plaquettes adja¬ 
cent to Z, l respectively. Similarly, 

Br=\(l+Y[a’ p ) , £ e =I(l + ]>!) ( 10 ) 

pec pec 


Bl = \(l + Bl)-P c , B\= l -{l+Bl)-P & . (12) 

Here P c is a projector that projects onto states obeying 
the closed membrane constraint in the neighborhood of 
the cube c. More specifically, 

p c = U A i U A i ( 13 ) 

l ^ c i± c 

where the first product runs over the twelve edges l of 
c and the second product runs over the six edges l that 
are perpendicular to the six plaquettes of the cube c. 
Similarly, 

p* = Il A i U A ' ( 14 ) 

lec l±a 

The Bl and B\ operators have a similar structure as m, 
but their matrix elements have some additional phase 
factors: 


Bl\X b ,X r ) = (-l) m 'i n '\X b + c,X r ), 
B\\X b ,X r ) = (-l) m H n *\X b ,X r + c). (15) 
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Here ra c , mg, n c , ng are integer-valued functions of 
(Xb,X r ), which we will now define. In fact, we will only 
define m c ,n c for membrane states (Xb,X r ) that satisfy 
P c = 1, and we will only define metric for states with 
Pc =3 1 ; it suffices to discuss these subsets of states since 
the projectors P c ,Pc ensure that the matrix elements of 
B\ and B\ vanish for all other states. 

We begin with m c . To define the value of m c for some 
(Xfr, X r ), consider the intersection between the set of red 
membranes X r and the cube c . This intersection defines 
a collection of closed loops that live on the surface of the 
cube c (Fig. 0^). (The fact that the loops are closed 
follows from the closed membrane constraint P c = 1 ). 
We will call these loops “red loops.” The integer m c is 
defined to be the number of these red loops, i.e., 

m c (X^, X r ) = #{red loops on c} (16) 

(See Fig. 0]I for an example). 

To define the value of n c for some (Xb,X r ), consider 
the intersection between the set of blue membranes X 5 
and the cube c. Because the blue membranes X 5 and 
the cube c both live on the same cubic lattice, the in¬ 
tersection X\) n c consists of 2D regions instead of loops. 
These 2D regions live on the surface of the cube c (Fig. 
0 )d). Consider the boundary of these 2D regions. This 
boundary consists of a collection of closed loops. We will 
call these loops “blue loops.” The integer n c is defined 
to be the number of intersections between the blue loops 
and the red loops defined above: 

n c (Xb, X r ) = ^{blue-red intersections on c} (17) 

(See Fig. 0]I for an example). Notice that n c is always 
even. 

The integers rric and ng are defined in an identical way 
but with the colors reversed: 

mc(I&, X r ) = #{blue loops on c} (18) 

and 

nc(Xi), X r ) = ^{red-blue intersections on c} (19) 

Here, the blue loops are defined by the intersection be¬ 
tween the blue membranes X 5 and the red cube c, while 
the red loops are the boundaries of the intersections be¬ 
tween the red membranes X r and the cube c. 

Alternatively, we can express m c , n c , ng in the spin 
language. In this language, the operator m c can be writ¬ 
ten as a function of the twelve spins cr| x ,..., <rj i2 that are 
closest to c (Fig. [5^), and similarly for rag: 

m c = 

mc = f{e z pi ,a* 2 ,...,(r z pi2 ). (20) 

Here / is a finite polynomial, but the explicit form of / is 
not illuminating so we do not show it here. Likewise, we 
can write n c as a function of the 6 blue spins cr^ ,..., cr | 6 



FIG. 5. (a) The integer m c is a function of the 12 spins cr| 
that are closest to the cube c. (b) The integer n c is a function 
of the 6 blue spins and 12 red spins cr| around the cube c. 

and 12 red spins, ,..., <rj i2 surrounding the cube c (Fig. 
Eb), and similarly for n 5 : 

n c = g(W z p },Wf,}), 

n e ~ g(Wp},Wp})- ( 21 ) 

As above, the expression for g is complicated so we do 
not show it here. 

While the Hamiltonian Hq is essentially identical to 
the well-known 3D toric code mode l 21 1 23 ~ —, some read¬ 
ers may be curious about the origin of the Hamilto¬ 
nian Hi. As we mentioned previously, this Hamiltonian 
has been designed to have a particular ground state - 
namely, | T 1 ). The state | T 1 ) is in turn motivated by the 
ground state of the spin model in Ref. [27| which describes 
a nontrivial symmetry-protected topological phase with 
Z 2 x Z 2 symmetry. Thus, the main question is how to 
design an exactly soluble Hamiltonian with a particular 
ground state, |Ti). One way to do this is to follow a 
similar approach to the string-net construction of Ref. 
I23I . We recall that the string-net models of Ref. [H can 
written as a sum of two types of operators — a Q op¬ 
erator and a B operator. The Q operator prefers closed 
string configurations and the B operator creates a closed 
loop. While the form of the Q operator is simple and in¬ 
tuitive, the B operator is more complicated, and its ma¬ 
trix elements are obtained by fusing a closed loop onto 
the lattice using 2D local rules. The Hamiltonian Hi 
can be constructed using a similar approach: in this case 
the Hamiltonian is built out of A operators which prefer 
closed membrane configurations and B operators which 
create a cube. Similarly to the string-net models, the pre¬ 
cise form of the B operator can be obtained by writing 
down local rules obeyed by |\Eb), and then fusing a cube 
onto the lattice using 3D local rules. This is one way to 
obtain Hi. Another way is to use membrane operators: 
as we will see in section UlI C 41 the B c , B g are examples 
of spherical membrane operators (|35I38|) . 

E. Properties of the Hamiltonians 

The two Hamiltonians Ho, Hi have many nice prop¬ 
erties. One property is that all the operators 
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[ At , A - t , Bj?, B?. B ], B \} are Hermitian so that Ho and 
Hi are also Hermitian. Indeed, the Hermiticity of 
{Ai,Aj, B®, B®} is clear from the definitions ([9110]) . while 
the Hermiticity of B\ and B\ can be seen by noting that 
the matrix elements of B\ and B\, as defined in (fl5|) , are 
both symmetric and real. 

Another important property is that all the terms in 
Ho, H\ commute with one another, so that the two mod¬ 
els are exactly soluble. For the case of H 0 , the fact that 
{Ai, Aj, B®, all commute with one another follows 
easily from the definitions (19IT01) of these operators. For 
the case of Hi, more work is required to verify this com¬ 
mutativity: while simple algebra shows that 

[A l ,A l ,] = [A i ,A i ,] = 0 (22) 

and 

[A U B 1 C ] = [A h B l c ] = [A U B\] = [A h B\] = 0, 

(23) 

the fact that 

[Bl,Bl,] = [B\,B\,} = [BlBl,] = 0 (24) 

is not obvious. We leave the derivation of the latter iden¬ 
tity (1M1) to Appendix [AJ 

A third property of Ho, Hi is that the operators 
{ Ai, Aj, B®, F>?, B\, B\} have eigenvalues 0 or 1: 

aha t ,b%b%b\,b\= 0,1. (25) 

(The first two eigenvalue spectra can be derived from 
Af — Ai, etc.). A final property of these models is that 
|To) is an simultaneous eigenstate of {Ai, Af, B®, B®} 
with 

ai = a i = \P C = 6? = 1. (26) 

Similarly, |Ti) is a simultaneous eigenstate of 
{Ai,A~ v Bl,B\} with 

ai = a t = b\ = b\ = 1. (27) 

A derivation of the relations ([26]) and (j27j) is given in 
Appendix [A] 

F. Solving the models 

In this section, we show that Ho and Hi are gapped, 
and that |4/o) and |Ti) are ground states of these Hamil¬ 
tonians. We begin with Hq- To find the energy spectrum 
of Ho, recall that the operators {Ai, A- { , B®, commute 

with one another and can therefore be simultaneously di¬ 
agonalized. Let us label these simultaneous eigenstates 
by |a*, a^, 6°, 6?) where ai, a^b^b® denote the eigenval¬ 
ues. It is clear that these states are energy eigenstates 
with energy 

£ = (28) 

i i c c 


Now, since the eigenvalues {ai, a^b^b®} take values in 
0 , 1 , it follows that the ground state(s) of Ho have ai = 
a i = h® = b® = 1 , while the excited states have at least 
one equal to 0. We conclude that there is a 

finite energy gap, A = 1 , separating the ground state(s) 
and excited states. Furthermore, we can see that |To) 
is a ground state of Ho since it obeys (|26|) . All that 
remains is to determine the ground state degeneracy of 
Ho . This degeneracy depends on the global topology of 
the system on which Ho is defined. In an infinite non¬ 
periodic geometry, one can show that Ho has a unique 
ground state — namely, |To). On the other hand, in 
periodic (3D torus) geometry, it can be shown that the 
ground state degeneracy is 4 3 = 64. These degenerate 
ground states are characterized by different parities of 
noncontractible red and blue membranes along the x,y 
and z directions of the 3D torus. 

Now let us consider Hamiltonian H\. Similarly to Ho, 
we can simultaneously diagonalize {Ai, Aj, B\, B\} and 
label their eigenstates as \ai,ai,b\,b\). Following the 
same reasoning as above, we conclude that there is a 
finite energy gap, A = 1 , separating the ground state(s) 
and excited states. Furthermore, we can see that |Ti) is 
a ground state of Hi since it obeys (EH). Like Ho, it can 
be shown that Hi has a unique ground state in an infi¬ 
nite non-periodic geometry (|Ti)), and the ground state 
degeneracy is 64 in a 3D torus geometry. 

G. Particle-like and loop-like excitations 

If we examine equation (l28l) . we can see that Ho sup¬ 
ports both particle-like and loop-like excitations. An ex¬ 
ample of a particle-like excitation is a point defect where 

or 6 ? is equal to 0 instead of the ground state value of 
1 . Likewise, an example of a loop-like excitation is a line 
defect along which ai or cq is equal to 0 instead of the 
ground state value of 1. These line-like defects always 
form closed loops: to see this, note that the quantum 
numbers ai and obey the local constraints 

JJ(l-2a,) = l , IL 1 - 2 ^ 1 ( 29 ) 

l € 3 les 

for every site s in the cubic lattice and s in the dual cubic 
lattice. Here, the first product runs over the six edges l 
that are adjacent to 8 and similarly for the second prod¬ 
uct. The above constraints guarantee that each site s is 
adjacent to an even number of line-line defects, so that 
the defects always form closed loops. (To derive these 
constraints, note that they follow from the corresponding 
operator identities, rizes(l ~ 2 Aj) = — ^Aj) = 1 

which in turn follow from the definition ®). In exactly 
the same way, one can see that the Hamiltonian Hi sup¬ 
ports particle-like excitations with b\ or b\ is equal to 
0 , and loop-like excitations where a\ or is equal to 0 
along some closed loop. 

Below we will see that in both models, these particle- 
like and loop-like excitations have nontrivial braiding 
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statistics with one another and have similar properties 
to the charges and vortex loops in Z 2 x Z 2 gauge the¬ 
ories. In view of this connection, we will refer to the 
particle-like excitations with b® or b\ equal to 0 as “blue 
charges” and the excitations with b 9 or b\ equal to 0 as 
“red charges.” Similarly, we will refer to the loop-like 
excitations with a\ = 0 as “blue vortex loops” and the 
excitations with = 0 as “red vortex loops.” In the 
following, we will study these charge and vortex-loop ex¬ 
citations in more detail, with a focus on their topological 
properties. 


III. EXCITATIONS AND THE ASSOCIATED 
CREATION OPERATORS 

In this section, we find operators that create the charge 
and vortex loop excitations of Ho and Hi. These oper¬ 
ators are useful because their commutation algebra con¬ 
tains information about the braiding statistics of the as¬ 
sociated particles and loops. 


A. General picture for string and membrane 
operators 

In general, topologically nontrivial particle excitations 
cannot be created using local operators. Instead, the eas¬ 
iest way to create these excitations is to use string-like op¬ 
erators. In the following sections we will find string-like 
creation operators for each topologically distinct charge 
excitation a in Ho and H\. We will denote these opera¬ 
tors by W a {P ) where P is the path along which the string 
operator acts. These operators satisfy two key properties. 
First, if P is an open path, then when W a (P ) is applied 
to the ground state |\k), it creates an excited state |\£ ea .) 
with two charge excitations a at the two ends of P: 

W a (P) |*> = |*e*>. 

(Here the excited state \^/ e x) only depends on the end¬ 
points of P and not on the choice of path). Second, if P is 
a closed path, then W a (P) does not create any excitation 
at all: W a (P) |tf> oc |tf). 

In addition to these string operators, which create 
charge excitations, we will also find membrane operators 
that create vortex loop excitations. That is, for each 
topologically distinct vortex loop excitation a in Ho and 
Hi, we will find a corresponding membrane creation op¬ 
erator, which we will denote by M a (S ) with S being the 
surface where the membrane operator acts. These mem¬ 
brane operators satisfy similar properties to the string 
operators. First, if S' is a cylindrical surface, then when 
M a (S ) is applied to the ground state |\I/) it creates an 
excited state with two loop excitations at the ends of the 
cylinder: 

M a (S)\9) = \9 ex ). 


(Similarly to before, the excited state \^f ex ) only depends 
on the two boundaries of S and not on the choice of sur¬ 
face that joins them). Second, if S is a toroidal surface, 
then M a (S ) does not create any loop excitations at all: 
M a (S) |*>oc|tf). 

We can also consider membrane operators with other 
topologies beyond the cylinder and torus case: for exam¬ 
ple, later we will discuss spherical membrane operators. 
However, the structure of membrane operators is differ¬ 
ent for different topologies^ and thus each case has to 
be treated separately. Here we will focus on cylindrical 
and toroidal cases as they are sufficient for our purposes. 

The string and membrane operators have simple phys¬ 
ical interpretations: the string operator W a (P) describes 
a process in which a pair of charge excitations is created 
and then moved to the two ends of the path P. Like¬ 
wise, the membrane operator M a (S ) describes a process 
in which a pair of loops is created and then moved to the 
two ends of the cylinder S. On the other hand, if P is 
a closed path then W a (P ) describes a process in which 
a pair of charge excitations is created and then moved 
around P and annihilated with each other. Similarly, if 
S' is a toroidal surface then M a (S) describes a process 
in which a loop-antiloop pair is created, moved around S 
and annihilated with one another. 


B. String and membrane operators for Ho 

1 . String operators 

It is easy to find string operators W® , W® that cre¬ 
ate the blue and red charge excitations of Ho. These 
operators are given by 

WS(P) = n W r(P) =U a l ( 3 °) 

p_LP p±P 

where P is a path on the dual cubic lattice in the first 
expression, and a path in the cubic lattice in the second 
expression. The two products run over plaquettes that 
are perpendicular to these two paths. Equivalently, in 
the membrane representation, W 6 ° and W® are given by: 

W 6 °(P) \X b ,X r ) = (-l) N *\X b ,X r ), 

W?(P) \X b ,X r ) = (~l) Nr \X b ,X r ) (31) 

where N & and N r are the number of blue and red mem¬ 
branes that cross the path P. 

Let us now verify that when P is an open path, W®(P) 
creates blue charge excitations at the two endpoints of P. 
From the definitions it is easy to see that W®(P) 
commutes with all the terms in the Hamiltonian Ho ex¬ 
cept for and Sj? 2 where ci and C 2 are the two cubes 
at the endpoints of P. These two operators anticommute 
with W 6 °(P) rather than commute. We conclude that 
W®(P) |^o) is a simultaneous eigenstate of Ai,A$, B®, P? 
with eigenvalues 6^ = b® 2 = 0 and all other eigenvalues 
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equal to 1. Hence, W®(P) |T 0 ) contains two blue charge 
excitations at the endpoints of P. A similar argument 
shows that W® also creates red charge excitations at the 
endpoints of P. 


2. Membrane operators 

It is also easy to find membrane operators M^, M® 
that create the blue and red vortex excitations of Hq. 
These operators are given by 

<(S) = n o%, M° r (S ) = n a? (32) 
pes pes 

where S is a surface made up of plaquettes p living in 
the cubic lattice in the first expression, and a surface 
consisting of plaquettes p in the dual cubic lattice in the 
second expression. Equivalently, we can express and 
in the membrane representation as: 

Mg(S)\X b ,X r ) = | X b + S,X r ), 

M?(S)\X b , X r ) = \X b , X r + S) (33) 

where X b + S denotes the Z 2 addition operation defined 
below equation (j8j). 

We now check that when S is a cylindrical surface, 
Mff(S) creates blue loop excitations at the two ends of 
S. To establish this fact, we note that M®(S) commutes 
with all the terms in the Hamiltonian Ho except for 
when l lies along the two boundaries of S. These srfi 
operators anticommute with M®(S) rather than com¬ 
mute. It then follows that M®(S) |To) is an eigenstate of 
Ai, Af, 13®, Bp with eigenvalues ai = 0 along the bound¬ 
aries of S, and all other eigenvalues equal to 1. This 
establishes the claim. The same argument applies to 
Af°(S). 

It should be noted that the excitations created by M® 
and M r ° are not the most general possible vortex loop 
excitations. In fact, the Hq model supports three other 
types of blue vortex loops and three other types of red 
vortex loops. These excitations can be obtained by at¬ 
taching either a red charge, a blue charge, or both a red 
charge and a blue charge, to the vortex loops created by 
M® and M ®. We will label the four types of blue vortex 
loops by (6, qb,q r ), and the four types of red vortex loops 
by (r, qb,q r ) where q b , q r take values in 0,1. In this label¬ 
ing scheme, the excitations created by M® and M® are 
denoted by (6,0,0) and (r, 0,0), with the other excita¬ 
tions labeled according to the amount of charge attached 
to them. As we will see later, these loop excitations are 
all topologically distinct from one another in the sense 
that they have different braiding statistics. 

It is easy to find membrane operators that create these 
more general vortex loop excitations. For example, con¬ 
sider the operator M®(S) in the case where S is a cylin¬ 
der. When this operator is applied to the ground state 
| T 0 ), it creates two vortex loops of type (6, 0, 0) at the two 


ends of the cylinder. We can modify this operator to cre¬ 
ate other vortex loops (6, qb,q r ) by multiplying M®(S) by 
the string operators W 6 °(P) or W?(P') or Wg(P)-W?{P') 
where the paths P, P' run along the length of the cylin¬ 
der. In the same way, we can construct membrane cre¬ 
ation operators for the red vortex loops (r, qb,q r ) by mul¬ 
tiplying M®(S) by appropriate string operators. 

Before concluding, we should mention that there is yet 
another type of vortex loop excitation which we have not 
discussed, namely a composite of a red and blue vortex 
loop. This type of loop excitation can be obtained by 
fusing together a red and blue loop, and again it comes 
in four subtypes which we can denote by (rb, qb,q r )- I n 
what follows, we will generally ignore this additional kind 
of vortex loop excitation since its braiding statistics prop¬ 
erties are completely determined by the properties of the 
individual red and blue loops. 


C. String and membrane operators for H 1 

In this section, we construct string and membrane op¬ 
erators for the Hamiltonian Hi. The string operators are 
relatively easy: following the same analysis as above, it is 
simple to check that the same string operators that cre¬ 
ate the blue and red charge excitations in the Ho model 
can also be used to create the charge excitations in the 
Hi model. That is, just as in (J3T|) , we can set 

W b \P)\X b ,X r ) = (-l) N >\X b ,X r ), 

WpP)\X b ,X r ) =(-l) N r\X b ,X r ) (34) 

where Nb and N r are the number of blue and red mem¬ 
branes that cross the path P. 

We will need to do more work to construct membrane 
operators for Hi. The rest of this section is devoted to 
this problem. 


1. Cylindrical membrane operators for blue vortex loops 

We begin by finding a membrane operator M^(S) that 
creates blue vortex loops. To proceed, it is convenient 
to work in the membrane basis. Let S be a cylindrical 
surface made up of plaquettes p living in the cubic lat¬ 
tice. We need to define the action of M^(S) on a general 
membrane state \Xb,X r ). A natural guess, inspired by 
©, is that M£(S) should act as 

M£(S)\X h ,X r ) = f b (X b ,X r ,S)\X b + S,X r ) (35) 

for some complex-valued function fb{X b ,X r , S). Indeed, 
it is easy to see that any operator of this type has the 
property that it anticommutes with the srf\ terms that 
he along the two boundaries of S and it commutes with 
every other and ^ term. Hence, any operator of this 
type has the property that if we apply it to the ground 
state |Ti), it will create blue vortex loop excitations at 
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FIG. 6. (a) A blue cylinder S intersects with a red membrane 
and two blue membranes (only part of S is shown for clarity). 
One of the blue membranes is incident from above S and the 
other from below, (b) We represent the intersections as a 
picture drawn on the cylinder S. In this case, the picture 
consists of a red loop, and two blue regions — one with a 
solid boundary and one with a dotted boundary. 


the boundaries of S. The problem is that most of these 
operators also create many other excitations along the 
surface S since most of these operators do not commute 
with B\ and B\. Thus, our task is to choose fb appro¬ 
priately so that M£(S) does not create any excitations 
except at the boundaries of S. 

We now describe one choice of fb that does the job. Be¬ 
low, we will just present the definition of fb without any 
motivation. Later, in section UlI C 51 we will explain why 
this choice works and we will provide some motivation as 
to where it comes from. 

First, we set fb(Xb,X r , S) = 0 if the membrane con¬ 
figuration (Xb,X r ) violates the closed membrane con¬ 
straint anywhere in the neighborhood of S. More pre¬ 
cisely, fb = 0 if any link /, l that touches or intersects S 
is adjacent to an odd number of occupied plaquettes. On 
the other hand, if every such link is adjacent to an even 
number of occupied plaquettes, then fb(Xb,X r , S ) is de¬ 
fined in terms of the two intersections Xb D S and X r D S. 
The dependence of fb on Xb D S and X r D S is compli¬ 
cated, and in order to explain it, we first describe a way 
to represent Xb D S and X r C\S as a picture drawn on the 
surface of the cylinder S. The picture we will draw con¬ 
sists of a collection of red lines and blue regions with the 
red lines denoting the places where membranes in X r in¬ 
tersect 5, and the blue regions denoting the places where 
membranes in Xb intersect/overlap with S. We note that 
the boundaries of the blue regions correspond to places 
where a membrane in X 5 is incident upon the surface S. 
If this membrane is incident upon S from above, we will 
draw the corresponding boundary as a solid line, while if 
it is incident from below, we will draw the boundary as a 
dotted line (Fig. [ 6 ]). (Note that for this step to be well- 
defined, one needs to specify a convention for what side 
of S is defined as “above” S and what side is defined as 
“below” S. Equivalently, one needs to choose a normal 
vector to S). Putting this all together, the intersections 
XbC\S and X r n S can be represented by a picture, drawn 
on the cylinder S', of the form shown in Fig. [7] 

The value of fb(Xb,X r , S) is completely determined 
by the corresponding picture. Thus, fb can be thought 
of as a complex-valued function defined on pictures. All 
that remains is to specify this function. For reasons that 


FIG. 7. A typical picture representing the intersections of an 
(unlinked) blue cylinder with red and blue membranes. Here 
we draw the cylinder as a rectangle with top and bottom 
identified and the left and right being the two ends of the 
cylinder. 


will become clear later on, it is most natural to describe 
fb implicitly, through local constraint equations, rather 
than providing an explicit formula. More specifically, we 
define fb in terms of the constraint equations 

fb = -fo (■)» ( 36a ) 

h Q) (Q = -fb > ( 36b ) 

h (O) = fh ( ' } ’ (36c) 

fb () = h , (36d) 

fb (V "}) = fb (•), (36e) 

h () () h • (36f) 

fb Qs) =/b 0> O’ (36g) 

fb 0$) = 0 0 ’ (36h) 

fb ^ = - fb ^ , (36i) 

h (X ) = ~ h (X) m 

where it is understood that the value of fb depends only 
on the topology of the picture. That is, any two pictures 
that can be smoothly deformed into one another have the 
same value of fb . The meaning of the above constraint 
equations is as follows. The first equation (I36ap states 
that the value of fb for a picture with a closed red loop is 
equal to the value of fb for the same picture, without the 
loop, up to a factor of —1. Equation (I36b[) states that two 
pictures that differ by the recoupling of red curves have 
values of fb that differ by a factor of —1. The remaining 
equations have a similar meaning. 

In the above constraint equations, we have left out the 
shading of the blue regions. But it should be understood 
that the blue regions are shaded in a consistent way on 
both sides of the equations. For example, (|36d [) repre- 
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TABLE I. The function fb that defines the (unlinked) blue 
cylinder operator is completely determined by its values on 
four basic pictures that are drawn on the surface of the cylin¬ 
der. Here the blue cylinder is represented by a rectangle with 
upper and lower edges identified and with the left and right 
being the two ends of the cylinder. The corresponding values 
of fb are shown below the pictures. The integers q r ,qb = 0,1 
define four different functions fb and therefore four different 
cylinder operators. 


sents two constraints with different shadings 

A ()() =A (X) > 

TONX)- 

Similarly, Q36jD also represents two constraints with dif¬ 
ferent shadings 



To complete the definition, we impose a boundary con¬ 
dition on fb which states that fb = 0 if any of the red 
lines or blue regions touch the boundaries of S. With this 
boundary condition and the above constraint equations, 
the value of fb on any picture can be related to one of 
the four “basic” pictures shown in Table [fl Thus, once 
we specify the value of fb on these basic pictures, we will 
have completely specified fb and therefore (S'). 

The values for fb on the basic pictures are shown in the 
bottom row of Table HJ These values are parameterized 
by two integers qb,q r E {0,1}. Hence our construction 
actually defines four different membrane operators. In 
principle, we should label these operators by M^ b qb q ^ 
to make the dependence on qb,q r explicit. However, this 
notation is cumbersome so we will denote the membrane 
operators by with the understanding that MI is not 
fully defined until we specify qb,q r - Similarly to the Hq 
model, these four different blue membrane operators cre¬ 
ate four different types of blue vortex loops. These vortex 
loops differ from one another by the amount of charge 
that they carry and we will label them by (b,qb,q r )- 

At this point, we have fully defined the membrane op¬ 
erator M^(S). All that remains is to show that this mem¬ 
brane operator has the required properties. That is, we 
need to show that M^(S) creates blue vortex loops at its 
two boundaries and nothing else. We will establish this 
fact in section UlI C 51 


2. Cylindrical membrane operators for red vortex loops 

To construct cylindrical membrane operators that cre¬ 
ate red vortex loops, we follow exactly the same recipe 
as above but with the roles of “red” and “blue” reversed. 
First, we define 

M}{S)\X h ,X r ) = f r (X b ,X r ,S)\X b ,X r + S) (37) 

where f r {Xb,X r , S) is a complex-valued function. We 
then define f r {X b , X r , S ) in terms of the two intersecting 
sets Xb fl S and X r D S. As in the definition of fb , we 
represent Xb D S and X r D S in terms of a picture drawn 
on the surface of the cylinder S, and we think of f r as 
a function defined on such pictures. In this case, each 
picture consists of a collection of blue lines and red re¬ 
gions, where the blue lines denote the intersection X^nS, 
and the red regions denote the intersection X r D S. We 
define the value of f r on each picture using local con¬ 
straint equations which are identical to Eqs. ([36]) . but 
with the red and blue colors reversed. With these con¬ 
straint equations, the value of f r on any picture can be 
related to one of the four basic pictures shown in Table 
[T] (with the colors reversed). We then define the value of 
f r on these basic pictures, just as in Table U but with the 
colors reversed and qb and q r exchanged. This procedure 
completely specifies f r and therefore M*. Note that, like 
M^, this construction actually gives four different mem¬ 
brane operators that are parameterized by two integers 
q b ,q r E {0,1}. These membrane operators create four 
different types of red vortex loops, which we denote by 

(r, qb,q r )- 


3. Cylindrical membrane operators for linked loops 

We now have all the tools we need to create an ex¬ 
cited state with a pair of vortex loops: to do this we 
simply apply one of the above cylindrical operators to 
the ground state |Tq). But what if we want to build 
a state with more than two vortex loops? One might 
try to make such a state by applying multiple cylin¬ 
der operators to the ground state. This approach will 
work fine if the state we are trying to build does not 
contain any linked loops. However, it will fail if any of 
the loops are linked: the above membrane operators are 
simply incapable of creating excited states with linked 
loops. To see this, imagine we first apply a cylindrical 
membrane operator M^(S) to the ground state, obtain¬ 
ing M 6 1 (5 ')|\Eh). This state contains two blue vortex loop 
excitations located at the two boundaries of the cylin¬ 
der S. Now, suppose we apply another cylinder operator 
M b (S') where S' is linked with one of these loops. One 
might hope that the result would be a state with linked 
vortex loops. Unfortunately, however, Ml {S') simply 
annihilates this state: Ml {S')Ml {S) |\hi) =0. To see 
this, note that Ml{S) |\Li) is a superposition of mem¬ 
brane states |Xfc,X r ), where all the membranes in X b 
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FIG. 8. Two typical pictures representing the intersections 
between (linked) blue cylinders and red and blue membranes. 
Here we draw the cylinder as a rectangle with top and bottom 
identified. Panel (a) shows a typical picture for a blue cylinder 
linked with a red base loop. Panel (b) shows a typical picture 
for a blue cylinder linked with a blue base loop. 


and X r form closed surfaces, except at the two bound¬ 
aries of S', where a blue membrane terminates. Since 
S' is linked with one of the boundaries of S, it follows 
that each of the two boundaries of S' must intersect with 
blue membranes an odd number of times. In particular, 
each boundary must have at least one intersection with 
a blue membrane. But if we look back at the boundary 
conditions for /&, we see that any state of this kind has 
fb{Xb, X r , S') = 0 , so any state of this kind is annihilated 
by Ml {S'). ' 

We therefore need to find another set of membrane 
operators to create linked loops. We now define such 
membrane operators. Before doing this, we first need to 
explain the basic idea behind these operators. In general, 
we will be interested in excited states containing a collec¬ 
tion of vortex loops that are all linked with a single “base 
loop” which can be either blue or red. We will build such 
states with the help of two different cylindrical membrane 
operators — one for each type of base loop. The operator 
associated with the blue base is designed to be applied to 
states containing a blue loop that links with the cylinder 
S: when the operator is applied to such a state, it creates 
two vortex loops at the ends of the cylinder S which are 
linked with the original blue loop. On the other hand, 
when it is applied to states that do not contain a blue 
loop that links with the cylinder 5, it simply annihilates 
them. The operator associated with the red base has 
a similar property. This structure is reminiscent of the 
cylinder operators constructed in the previous two sec¬ 
tions. In fact, the previously constructed operators can 
be thought of as cylindrical membrane operators associ¬ 
ated with a trivial base: they only create loop excitations 
when applied to states without any loops linked with S. 

With this picture in mind, we now construct cylindrical 
membrane operators for a red base loop and a blue base 
loop. We begin by defining a cylindrical membrane op¬ 
erator that creates blue vortex loops that are linked to a 
red base loop. This operator is defined in almost exactly 
the same way as the blue membrane operator M£(S) de¬ 
fined above. The only difference is that we change the 
boundary condition on fb so that it only takes nonzero 
values on pictures in which a single red line touches the 
boundary of 5 at a fixed position yo , and no blue regions 


eft 73 

V ► V 


‘ - 1 

A 

C 



3 : 


fb 

1 

i • e 

iTrq r 

i ■ e inqb 

—e 

,in{qr+q b ) 

d A/ 




h 


5 


fb 

1 

e i7rq r 

e ^q b 

e iTr(q r +q b ) 


TABLE II. The function f b that defines the linked blue cylin¬ 
der operator is completely determined by its values on four 
basic pictures that are drawn on the surface of the cylinder. 
The top panel shows the four basic pictures for a blue cylinder 
linked to a red base loop, while the bottom panel shows the 
pictures for a cylinder linked to a blue base loop. The cor¬ 
responding values of fb are shown below the pictures. Here 
q r ,qb — 0,1 define four different functions fb and therefore 
four different cylinder operators for each panel. 


touch the boundary. (See Fig. [H^i for a typical picture of 
this type). Given this boundary condition and the con¬ 
straint equations ([36]) we can reduce any such picture to 
one of the four basic pictures shown in the top panel of 
Table [XT] We can therefore completely specify the mem¬ 
brane operator once we specify the value of fb on these 
four basic pictures. The values we choose are shown in 

Table HI 

We next describe a cylindrical membrane operator 
which creates a blue vortex loop which is linked to a blue 
base loop. Here, we choose a boundary condition on fb 
so that it only takes nonzero values on pictures in which 
a single thin blue region touches the boundary of S at a 
fixed position yo, and no red lines touch the boundary. 
(See Fig. [ 8)3 for a typical picture of this type). Given this 
boundary condition and the constraint equations <155)1 , we 
can reduce any such picture to one of the four pictures 
shown in the bottom panel of Table [Til The value of fb 
on these four pictures is shown in Table [III 

To complete the discussion, we need to explain how 
to construct cylindrical membrane operators that create 
red vortex loops linked to a blue or red base loop. These 
operators are defined exactly like the two operators de¬ 
scribed above, but with the roles of “red” and “blue” 
reversed. 

Before concluding, we make one comment about our 
notation: we will use the same symbol M£(S) to denote 
all of the blue membrane operators whether they are as¬ 
sociated with a red base (top panel of Table ITT1) , a blue 
base (bottom panel of Table HI]) . or no base at all (Ta¬ 
ble [IJ ). Similarly, we will denote all the red membrane 
operators by M}(S) independent of what base they are 
associated with. This abuse of notation will not intro¬ 
duce confusion since it will always be clear from context 
which operator we have in mind. 
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4- Toroidal and spherical membrane operators 

So far we have only discussed cylindrical membrane 
operators. We now discuss how to construct membrane 
operators with other topologies — in particular, toroidal 
and spherical operators. We begin with the toroidal op¬ 
erators. Roughly speaking, these operators can be ob¬ 
tained by connecting the two ends of our cylinder opera¬ 
tors. More precisely, we define toroidal operators exactly 
like the cylindrical operators, with only one modification: 
we replace the cylindrical boundary conditions on f b and 
f r with periodic boundary conditions in both directions 
of the torus. 

Like cylindrical operators, we can build different torus 
operators for different base loops. The blue torus oper¬ 
ator for a trivial base loop is defined by the values of 
fo shown in Table U while the blue torus operator for a 
red base loop is defined by the values shown in the top 
panel of Table [TTJ Finally, the blue torus operator for a 
blue base loop is defined by the values of /*, shown in the 
bottom panel of Table [TTJ Red torus operators for various 
base loops can be defined in the same way but with the 
red and blue colors exchanged. 

It is important to keep in mind that the two coordi¬ 
nates that describe the torus are not equivalent. One 
coordinate parameterizes the movement of the loop in 
space, while the other coordinate parameterizes the loop 
itself. Thus to define a torus operator, we not only have 
to specify the torus S', but we also have to specify which 
coordinate has which meaning. 

Let us now discuss spherical membrane operators. We 
can define blue spherical operators following the same 
recipe as the blue cylindrical operators. The only differ¬ 
ence is that in the spherical topology, every picture can be 
reduced to the vacuum or (empty) picture by application 
of the constraints ([36]) — in contrast with the cylindrical 
case where every picture can be reduced to one of the 
four pictures shown in Table U or Table m Thus, in the 
spherical case it suffices to define the value of fo on the 
vacuum picture. Here we define /^(vacuum) = 1 . Note 
that there is only one type of blue spherical operator — 
unlike the cylindrical or toroidal case where there are four 
types of operators parameterized by q r , G {0,1}. 

In fact, we already encountered spherical membrane 
operators in the definition of the Hamiltonian H\\ the 
operators B\ and B\ (H3 can be thought of small spheri¬ 
cal membrane operators associated with a unit blue cube 
S = c or unit red cube S = c. To see that Bl is a spher¬ 
ical membrane operator, note that the picture drawn on 
the surface of the blue unit cube has the property that 
the boundaries of the blue regions are always solid lines 
—that is, there are no boundaries that are dotted lines. 
The reason this is so is that all membranes in X 5 are inci¬ 
dent on the cube from the outside, rather than the inside. 
Now, since the picture on the surface of the cube does not 
contain any dotted blue lines, the local constraints ([36 [) 
simplify considerably: in fact, we can throw out equa¬ 
tions (I36e|) - ( |36gP and (|36il|36jD . The resulting equations 


can be solved explicitly, leading to the following formula 
for f b : 

f b (X b , X r , S ) = (-l)#f red 1o °p s 011 c > 

^#{blue-red intersections on c} (38) 

We can see that this formula, together with the definition 
([35]) agrees exactly with the definition of B\. 


5. Showing that the membrane operators have the required 
properties 

Having defined the membrane operators M£(S) and 
M^(5), we now show that they have the required prop¬ 
erties. We will focus on one case, namely blue membrane 
operators M£(S) that create unlinked loop excitations. 
The arguments for the linked case are similar. 

To begin, we show that if S is a torus then M^(S) does 
not create any excitations at all. That is, we show 

M£(S)\*i) (X (39) 

To establish this result, we first rewrite dm in a more 
convenient form. Multiplying both sides by (X b , X r | 
gives: 

(X b ,X r |M b 1 (5)|«' 1 ) oc <X 6 ,X P |tf 1 ). (40) 

Next, using the definition of M£(S) ([35]) . we can rewrite 
this as 

ft(X b ,X r ,S) • ^(W, + S,X r ) oc Tx(X 6 ,X r ). (41) 

To proceed further, we observe that Eq. m is equivalent 
to the relation 

/ 6 *(x 6 ,x r ,g).* 1 pr 6 + s,x r ) = *i(WW) ( , 2) 
r b {X' b ,X' r ,S)-y l {X' b + S,X' r ) *i{X' b ,X' r ) 1 ’ 

for any two closed membrane configurations Xb,X r and 
X' h ,X' r . Finally, we use the fact that f b is a pure phase 
to rewrite this equation as 

h{x' b ,x' r ,s) = 4/i(WW)*i(^ + ff,V) 
f b (x b ,x r ,s ) (xi„xoMx b + s,Xry [ ’ 

Our task is now to prove equation (1151) . First, we claim 
that it suffices to prove m for the case where (X^X' r ) 
and (X b , X r ) only differ locally — i.e. only differ in some 
small region. The reason it is enough to consider this 
case is that we can get from any membrane configuration 
(X b ,X r ) to any other configuration (X^X!^) by a series 
of local changes: 

(X b ,X r ) -► {X bl ,X rl ) -► (X b2 ,X r2 )---^ (Xix' r ). 

If (@3|) holds for each pair (X bi ,X H ),(^ 6 (i+i),V(i+i)) 
then once we multiply these relations together, we see 
that it automatically holds for (X^X!^) and (X b ,X r ). 
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FIG. 9. A blue toms operator Mb(S) acts on two slightly 
different membrane configurations (a) (Xb,X r ) and (b) 
(X(,jXr). Here, Xb = X' h is a blue sphere, while X r con¬ 
sists of two red spheres and X' r consists of a bigger red sphere 
obtained by merging the two red spheres in X r . 


There are two cases to consider: the region where 
(X b , X' r ) and (X bl X r ) differ may overlap S or not overlap 
S. In the second case, it is easy to see that (Iff3l) holds. 
Indeed, in this case the left hand side of (|43l) is equal 
to 1 since fb only depends on the membrane configura¬ 
tion in the neighborhood of S and hence fb(X b , X ' r , S) = 
fb{Xb , X r ,S). Similarly, the right hand side of (|43|) is also 
equal to 1 since the ratio Ti(X^,X')/Ti (Xb,X r ) only 
depends on the membrane configuration in the region 
where ( X' h ,X' r ) and (X b ,X r ) differ from one another. 

All that remains is the case where (X' h ,X' r ) and 
(. Xb,X r ) differ in a small region that overlaps S. In this 
case, the ratio on the left hand side of (|43|) is directly 
determined by the local constraint equations ([36]) for f b . 
The key point is that these constraint equations were 
chosen specifically so that a solution to these equations 
will automatically obey (li51) . This is easiest to see by 
example. Consider a blue torus operator M b (S) acting 
on the two configurations (Xb,X r ) and (X' h ,X' r ) shown 
in Fig. [9] Here X b = X' b consists of a single blue sphere, 
while X r consists of two red spheres and X' r consists of 
single red sphere obtained by merging together the two 
spheres in X r . In this case 


so that 


h{X' b ,X' r ,S) 

f b (X b ,X r ,S) 



= -1 


(44) 


according to (I36bll . On the other hand, it is easy to see 
that 


*l(WW) = *1 {Xh+^Xr) 

Vi(K x t) ’ 4-1 (X' b + S,Xl) 

using the explicit formula for 4/1 (f2ll . We conclude that 

(*£ + £,*;) = 

(x',x;)^!(x 6 + 5,x r ) 1 ' 

so that equation (1151) is satisfied in this case. 


The above example serves two purposes. First it 
demonstrates, in at least one case, that fb obeys equa¬ 
tion m- Second, it reveals where the constraint equa¬ 
tion (l36bl) comes from: it should be clear that we chose 
the factor of —1 in this equation specifically to ensure 
that (1151) was satisfied. Similarly, one can check that the 
other constraint equations ([36]) ensure that fb obeys (l43|) 
in other cases. 

This concludes our argument showing that the torus 
operator M b (S) does not create any excitations when 
applied to the ground state. Next, we need to show that 
the cylinder operator M b (S) does not create any excita¬ 
tions except at the two boundaries of S. This claim can 
be established using a similar argument to the torus case, 
but we will not repeat the derivation here. Instead we 
simply observe that the cylinder and torus operators look 
identical except near the two boundaries of the cylinder. 
Therefore it is intuitively clear that since the torus oper¬ 
ator does not create any excitations, the same must be 
true of the cylinder operator away from its boundaries. 

At this point, we have argued that the cylinder opera¬ 
tor M b (S) creates blue loop excitations at its two bound¬ 
aries and no other excitations anywhere else. However, 
we are not quite finished: we still need to verify one more 
property of the cylinder operator M b (S). Recall that Hi 
supports four topologically distinct types of blue loop 
excitations (6, qb,q r ) which differ from one another by 
attaching red and blue charges (see section IIIIB2D . We 
need to check that M b (S) creates exactly one of these 
excitations, and not a linear superposition of different 
types of excitations. In other words, we need to check 
that M b (S) creates loop excitations that are eigenstates 
of braiding measurements. 

To see that M^ 1 creates braiding eigenstates, we make 
use of a result from appendix [B] In that appendix, we 
show that the cylinder operators M b (S) are guaranteed 
to create braiding measurement eigenstates provided that 
fb is multiplicative in the sense that 

f b (X b ,X r ,SU S’)*= f b (X b ,X r ,S) ■ f b (X h ,X r ,S') (46) 

for any two cylinders S and S' that share a com¬ 
mon boundary and any membrane configuration (X*,, X r ) 
whose intersection with the common boundary obeys the 
appropriate cylinder operator boundary condition. 

In view of this result, we only have to show that fb 
obeys condition (0SJ). We go through this calculation in 
appendix [Bj and we show that fb does indeed obey (|46|) 
provided that the values of fb on the four basic pictures 
are those shown in Table H In fact, this is why we picked 
the particular values shown in that table: we chose those 
values to ensure that fb obeys equation 

To summarize, we have shown that the cylinder op¬ 
erators M b (S) create blue loop excitations at their two 
boundaries and no other excitations anywhere else. We 
have also shown that these blue loop excitations are 
eigenstates of braiding measurements. This establishes 
that the cylinder operators have all the required proper¬ 
ties. 
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D. Labeling scheme for loop excitations 



In this section, we discuss some subtleties related to 
the labeling of loop excitations. As we have emphasized 
previously, both the Hq and Hi models support four dif¬ 
ferent kinds of blue loop excitations and four different 
kinds of red loop excitations. We label the former by 
(b,qb,q r ) and the latter by (r, qb,q r ) where qb,q r = 0,1. 
In practice, we first assign labels to the membrane op¬ 
erators; we then assign labels to the loop excitations ac¬ 
cording to which membrane operator creates them. 

Now, an important question is how we choose which 
loop excitations are labeled by (6,0,0) and (r, 0,0). We 
will call these excitations neutral loops. Once we decide 
which blue and red loops should called neutral, the la¬ 
beling of all the other excitations is naturally fixed: for 
example, we assign the label (6, qb,q r ) to the loop that is 
obtained by attaching qb blue charge and q r red charge 
to the neutral loop, (6, 0, 0). 

For the case of unlinked loops, there is a natural choice 
for which loop should be called neutral: the neutral loop 
is the unique loop that can be shrunk down to a point 
and annihilated by a local operator. Thus, for the case 
of unlinked loops there is a canonical labeling scheme. 
This is the labeling scheme we use here. Indeed, in the 
case of Hq, we assign (6,0,0) to the loop created by the 
operator and one can readily verify that this loop 

is the unique blue loop that can be annihilated locally. 
Similarly, in the case of Pi, we assign (6, 0,0) to the loop 
created by the membrane operator defined in Table |T] with 
qb = qr = 0 and one can also check that this loop can be 
annihilated locally. (One way to see this is to note that if 
we shrink the two ends of the cylinder to make a sphere, 
the operator defined by Table [T] reduces to the spherical 
membrane operator defined by ®ED). 

In contrast, for linked loops, there is no canonical way 
to define which loops are neutral. Therefore, we make 
this assignment arbitrarily. In other words, for each base 
loop, we arbitrarily choose one of the blue linked loops 
and one of the red linked loops, and we call them neu¬ 
tral. In the case of Pi, this arbitrary choice enters in 
how we parameterize the values in Table HD which define 
the membrane operators and hence define the loop exci¬ 
tations. As we show in appendix [Bj the values in Table 
mi are obtained by solving certain algebraic equations. 
These equations have four different solutions, which de¬ 
fine four different membrane operators. All four are on 
an equal footing, but we simply picked one and labeled 
it by qb = 0, q r = 0. The labeling of the other solutions 
is then completely fixed. We could replace qb —>> qb + 1 or 
q r — »■ q r + 1, for either the top or bottom panel of Table 
im and the resulting table would define an equally valid 
labeling convention. 

The arbitrariness of our labeling scheme has an im¬ 
portant consequence: because we picked the labels in an 
arbitrary way, there is no sense in which an unlinked loop 
of type (6, qb,q r ) is the “same” type as a linked loop of 
type (6, qb,q r )- More generally, we cannot sensibly com- 


— ► < - 

p, p 3 

FIG. 10. The exchange statistics of the charge excitations 
can be computed from the commutation algebra (l47l) of three 
string operators acting on three paths Pi, P2, P3 that share a 
common endpoint. 

pare loops that are linked with different base loops. Each 
base loop effectively defines its own universe of excita¬ 
tions. 


IV. BRAIDING STATISTICS OF EXCITATIONS 

In the previous section, we constructed operators that 
create charges and vortex loops for the Po,Pi mod¬ 
els. In this section, we will use these operators to com¬ 
pute the braiding statistics of these excitations. We con¬ 
sider four types of processes: (1) processes involving two 
charges, (2) processes involving a charge and a vortex 
loop, (3) processes involving two vortex loops, and finally 
(4) “three-loop processes” involving two vortex loops that 
are both linked with a third loop. We find that the two 
models have the same braiding statistics except for the 
last process; when we investigate the three-loop braid¬ 
ing process we find a distinction between the two models 
which implies that they belong to distinct topological 
phases. 

A. Braiding two charges 

To compute the statistics of the charge excitations, we 
use the “hopping operator algebra” derived in Ref. [H. 
This algebra relates the exchange statistics of particle 
excitations to the commutation properties of the string 
operators that create these particles. To see how this 
works, let us consider the blue charge excitations in the 
Hq model. According to the hopping operator algebra, 
the exchange statistics of these excitations can be read 
off from the phase factor in the commutation relation 

wg (Pi)wg (P 2 ) w 6 ° (P 3 ) I ^ 0 } = 

e ie -W^P 3 )WH(P 2 )W°(P 1 )\^ 0 ) (47) 

where Pi, P 2 , P 3 are three paths arranged in the geome¬ 
try of Fig. flOland |4>o) is the ground state of Hq. That is, 
if 6 = 0, then the blue charges are bosons, while if 0 = 7r, 
the blue charges are fermions. (Other values of 0 are not 
possible in 3D). If we examine the form of the string op¬ 
erators, we can see that 0 = 0 since the string operators 
IF 6 0 (Pi), IF 6 0 (P 2 ), W^(Ps) (l30l) all commute with one an¬ 
other. We conclude that the blue charges are bosons in 
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FIG. 11. (a) The statistical phase associated with braiding a 
charge around a loop can be computed from the commutation 
algebra (|48l) of a string operator acting along a path P, and 
a membrane operator acting along a cylinder S. (b) The 
phase associated with braiding two loops around one another 
can be computed from the commutation algebra (HU of two 
membrane operators, one acting on a cylinder S , and the 
other acting on a torus S '. 


the Hq model. Using identical reasoning, we can see that 
the blue and red charges are bosons in both models. 


B. Braiding a charge and a vortex loop 

Next, we compute the statistical Berry phase associ¬ 
ated with braiding a charge around a vortex loop. More 
specifically, let us consider the Hq model and imagine 
braiding a blue charge around a blue vortex loop. The 
statistical phase 0 for such a process can be read off from 
the commutation algebra of the corresponding string op¬ 
erator and membrane operator 

M$(S)WS(P )|* 0 > = e ie ■ Wj}(P)Mg(S) |tt 0 > (48) 

where P is a closed path and S is a cylindrical surface, 
arranged as in Fig. ITTfa . To see where this relation comes 
from, note that the operator M b (S) describes a process in 
which two vortex loops are created and moved to the ends 
of the cylinder S, while the operator W b (P) describes a 
three-step process in which two charges are created out 
of the vacuum, one of them is braided around the path 
P, and then the two annihilated with one another. With 
this interpretation, the right hand side of ([48]) describes 
a process in which two vortex loops are created and then 
a charge is braided around P, while the left hand side 
describes a process in which a charge is braided around 
P first, and then two vortex loops are created. Clearly 
the phase difference between these two processes is the 
statistical phase associated with braiding a charge around 
a vortex loop. 

Remembering the definition of M b (S) and IU 6 0 (P), 
we can see that these operators anticommute with each 
other. We conclude that braiding a blue charge around a 
blue vortex loop results in a statistical phase of 7r. Sim¬ 
ilarly, we can see that M®(S) and W®(P) anticommute 
with each other so braiding a red charge around a red 
vortex loop gives a phase of n. On the other hand, if 
one braids a red charge around a blue vortex loop or 
a blue charge around a red vortex loop, the statistical 
phase vanishes since these operators commute with one 
another. 


A similar calculation for the Hi model gives identical 
results. Thus the two models share the same statistics 
between charges and vortex loops. We note that these 
statistical phases agree with the Aharonov-Bohm phases 
associated with braiding a charge around a vortex loop 
in Z2 x Z2 gauge theory. This is not a coincidence: as 
we mentioned in the introduction, the two models can be 
viewed as two different types of Z2 x Z2 gauge theories. 


C. Braiding two vortex loops 

1. Ho 

We now consider, for the Hq model, the statistical 
phase associated with braiding a vortex loop around an¬ 
other vortex loop (Fig. [TJl). For example, let us consider 
braiding a blue vortex loop around another blue vortex 
loop. Similarly to Eq. ([48]) . the statistical phase 6 bb can 
be computed from the commutation relation 

M°{S')M°{S) |tt 0 > = e<> ■ MS(S)M$(S') |# 0 ) (49) 

where S' is a cylinder and S' is a torus, arranged as in 
Fig. [Tib . Examining the definition of the membrane 
operators we can see that they commute with one 
another so that 0 bb = 0. Likewise, we can see that there 
is no statistical phase associated with braiding a red loop 
around a red loop or a red loop around a blue loop since 
the corresponding membrane operators all commute with 
one another. 

The above results apply to the blue vortex loop and red 
vortex loops labeled by (6,0,0) and (r, 0,0). One might 
also wonder about the braiding statistics of more general 
vortex loops (6, qb,q r ) and (r, <#>, g r ). The braiding statis¬ 
tics of these more general vortex loops can be computed 
using the same approach as above. The only difference 
is that the membrane creation operators are slightly dif¬ 
ferent in this case: they are obtained by multiplying M b 
(or M®) by one of the string operators W b (P ), W^(P') 
or W b (P) • W®(P'), where the paths P, P' run along the 
length of the cylinder. Substituting these modified mem¬ 
brane operators into Eq. a simple calculation shows 
that the phase associated with braiding a blue vortex 
loop (6, qb,q r ) around another blue vortex loop (6, q' b ,q' r ) 
is given by 


Qbb = A<ib + q' b )- ( 5 °) 

Similarly, the phase associated with braiding a red vortex 
loop (r,qb,q r ) around another red vortex loop (r,q' b ,q' r ) 
is 


e°r = 7T (q r + q’ r ) (51) 

while the phase associated with braiding a red vortex 
loop (r, qb,q r ) around a blue vortex loop (6, q bl q' r ) is 

0rb = n{qb + q' r ). (52) 
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These expressions have a natural interpretation in terms 
of Aharonov-Bohm phases: the first term in each for¬ 
mula is the statistical phase associated with braiding the 
charge on the first vortex loop around the flux of the sec¬ 
ond loop while the second term is the phase associated 
with braiding the flux of the first vortex loop around the 
charge on the second loop. 


2. H x 

Similarly to the Ho case, the statistical phase 6 bb as¬ 
sociated with braiding a blue vortex loop around a blue 
vortex loop in the Hi model can be read off from the 
commutation algebra of the corresponding membrane op¬ 
erators: 

MZ(S')MZ(S) |#i) = e< b ■ Ml(S)Ml{S') |*i>. (53) 

However the membrane operators M b are more compli¬ 
cated than their counterparts in the Hq model, so it is 
more difficult to apply Eq. ([53]) . To deal with this issue, 
we now derive a simpler version of Eq. (El) which is more 
convenient for our purposes. The first step is to rewrite 
Eq. ([53]) as 

\X)(X\V 1 ) = 

X 

e< b Y, M b( S ) M b(S')\X)(X\^ 1 ). 

X 


where \X) = \X b ,X r ). Next, we use the fact that 
Ml(S')Ml(S)\X b ,X r ) oc \X b + S + S',X r ), which im¬ 
plies that all of the states M b (S')M b (S)\X) that appear 
in the above sum are linearly independent from one an¬ 
other. It follows that the above equality must hold for 
each term separately. Hence, we must have 

Ml{S’)Ml{S)\X) = e<» ■ Ml(S)Ml{S’)\X) (54) 


for every membrane state \X) that has a nonzero ampli¬ 
tude in the state |Ti). 

We will now use (1M1) to compute the statistics of the 
loops in the Hi model. To this end, we set \X) = |0, 0), 
the “no-membrane” state. We then define | Y) = 15,0), 
\Y') = |S",0>, and | Z) = |S + S",0>. With this notation, 
the left side of ([54]) can be computed as 

MpS’)MpS)\X) = MpS') • f b (X,S)\Y) 

mf b (Y,S')f b {X,S)\Z). (55) 

Similarly, the term on the right side is given by 

Ml{S)Ml{S’)\X) = Ml(S) • f b (X,S')\Y’) 

= f b (Y',S)f b (X,S')\Z). (56) 

Comparing these two expressions with (154]) . we derive 


<b = f b (Y,S')f b {X,S) 

MY',s)f b (x,s'y 


( 57 ) 


To complete the calculation, we compute the value of f b 
for each of these configurations 

Mx,s) = f b ( 

f b (X,S') = f b (^ 
f b (Y,S') = f b ( 
f b (Y’,S) = f b ( 

Substituting these values into E3, we find that the sta¬ 
tistical phase associated with braiding a loop (6, qb,qb) 
around another loop (b,q b ,q b ) is 

Q b b = n{q b + q' b ). (59) 



= 1 , 

= 1 , 

; e *7 rq b f 

- e inqb' 


(58) 


We note that this result is identical to the statistical 
phase 0 bb ([50]) in the Ho model. Similarly, it is easy 
to check that the phase associated with braiding two 
red loops around one another agrees with ([51 ]) while the 
phase associated with braiding a red loop around a blue 
loop agrees with ([52]) . Thus, the two models share the 
same two-loop braiding statistics. In fact the agreement 
between the two models is not surprising, since the two 
models are equivalent to Z 2 x Z 2 gauge theories, and in 
such systems the two-loop braiding statistics always has 
an Aharonov-Bohm form, as explained in Ref. [3. 


D. Three-loop braiding 

Finally, we discuss the three-loop braiding statistics 
in the two models. Specifically, we consider a braiding 
process in which a loop is braided around another loop, 
while both are linked to a third “base” loop (Fig. [TJd). 
Unlike the other processes we have considered until now, 
we will see that the two models can be distinguished by 
their three-loop statistics. 


1. Ho 

First we compute the statistics in the Ho model. As 
in the two-loop case, the three-loop braiding statistics 
can be read off from the commutation algebra for the 
membrane operators. For example, the statistical phase 
0 bb r associated with braiding a blue vortex loop around 
a blue vortex loop, while both are linked to a red base 
loop can be obtained from 

M°(S')M°(S)\y ex ) = ■ M b (S)M b (S')\^ex) (60) 

where S is a cylinder, S' is a torus, and \^ ex ) is an excited 
state with a red vortex loop that links with both S and S' 
(Fig. IT2 ]) . Remembering that the membrane operators 
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FIG. 12. The statistical phase associated with a three-loop 
braiding process can be computed from the commutation alge¬ 
bra of two membrane operators, one acting along the cylinder 
S, and the other acting along the toms S '. Here both S and 
S' are linked with a base loop which lies along the boundary 
of the disk D. 


M b all commute with each other, we deduce that 6 bb r = 
0. In the same way, we can see that 0 bb b = 0. 

In the above calculation we implicitly assumed loops 
of the form (6,0,0). If we instead consider a process in 
which a general vortex loop (6, q b ,q r ) is braided around 
another loop (b,q' b ,q' r ) with either a red base loop or a 
blue base loop, then, just as in the two-loop case ([50|). 
one finds 


@bb,r ~ @bb,b ~ ^(Qb + Qb)- ( 61 ) 

Similarly, if we braid a red loop (r, q b ,q r ) around another 
red loop (r,q' b ,q' r ) with either base, the statistical phase 
is 


°rr,r = °rr,b = + q' r ). (62) 

Finally, braiding a red loop (r, qb, q r ) around a blue loop 
(6, q b ,q' r ) with either base gives the phase 

e °rb,r = 6rb,b = ' K {qb + q'r)- (63) 

The agreement with the two-loop statistics U5DK2D 
is to be expected since the membrane operators 
M b (S),M b (S') obey the same commutation algebra in¬ 
dependent of whether they act on the ground state |\ho) 
(as in the two loop case) or an excited state \^ ex ) (as in 
the three loop case). 


2. Hi 

We begin by computing the statistical phase associated 
with braiding a blue vortex loop around another blue 
vortex loop, while both are linked with a red base loop. 
As in the Hq case, the statistical phase 0l b r is given by 

Ml(S')Mt(S )|* ex ) = • M 6 1 (S)M 6 1 (S')l*e*> (64) 

where S is a cylinder, S' is a torus, and \^ e x) is an excited 
state with a red vortex loop that links with both S and S' 
(Fig. [12]) . Following the same logic as in section HV C 2[ 
it is straightforward to deduce an alternative and more 
convenient form of equation ([64]) : 

Ml{S’)Ml{S)\X) = e ie ^Ml{S)Ml(S’)\X) (65) 


where \X) = \X b ,X r ) is any membrane state that has 
nonzero overlap with |4/ e;r ). To compute 0 bbr , we set 
\X) = |0 ,D) where D is a disk with a boundary along 
the red vortex loop. We then define | Y) = \S,D), \Y') = 
| S' ,D), and | Z) = \S + S' ,D). With this notation, the 
left hand side of ([65]) can be computed as 

Ml(S')Ml(S)\X) = f b (Y,S')f b (X,S)\Z) (66) 

while the expression on the right hand side can be written 
as 


Ml(S)Mt(S’)\X) = f b (Y',S)f b (X,S')\Z). (67) 


The three-loop statistics is then given by 

ieL f b (Y,S')f b (X,S) 

f b (Y',s)f b (x,s>y 


( 68 ) 


To complete the calculation, we need to compute the 
values of f b on each of these configurations. We have 


Mx,s) = f b ( 

f b (X,S') = f b ( 

MY,s’) = f b ( 

f b (Y',S) = f b ^ 

where the ‘ —’ sign in the last equation comes from using 
the local rules (136] ) to reduce the picture to one of the 
basic pictures in Table [TTJ Substituting these values into 
([68]) . we find that the statistical phase associated with 
braiding a loop (b,q b ,q r ) around another loop (b,q b ,q' r ) 
while both are linked to a red base loop is given by: 

e bb,r = n + n(q b + q' b ). (70) 

In a similar manner, it is easy to show that the phase 
associated with braiding a red loop (r, q b ,q r ) around an¬ 
other red loop (r,q' b ,q' r ), while both are linked to a red 
base is: 



= 1 , 

= 1 , 

= i • . 


(69) 


8lr,r = 7T (qr + q'r)- ( 71 ) 

As another example, let us compute the phase associ¬ 
ated with braiding a red loop around a blue loop, with a 
red base. The analogue of equation ([65]) in this case is 

Ml{S')Ml{S)\X) = e id ’' b - r ■ Ml{S)M}(S’)\X) (72) 

where \X) = \X bl X r ) is any state that has nonzero 
overlap with \^ ex ) and \^ e x) is an excited state with 
a red vortex loop that links with both S and S'. Let¬ 
ting \X) = |0,£>>, \Y) = | S,D), | Y') = \0,D + S'), and 
| Z) = |S', D + S'), we derive 

i0lbr fr^SVkjXS) 

f b {Y\S)f r {X,S’Y 


(73) 
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Proceeding as before, we find 


h(x,s) 

fr(X,S') 

fr(Y,S') 

fb(Y',S) 


fr 

fr 

fb 




J 

V 

1\ 


1 , 

1 , 

e i ’ Kq ' b , 


i-e i7Tqr . 


(74) 


Substituting these values into (1731) , we conclude that the 
statistical phase associated with braiding a loop (r, qb,q r ) 
around another loop (b, q' b , q' r ) while both are linked to a 
red base loop is 


0r6 ) r~-|+ 7r (®>+<?r)- (75) 

So far we have only discussed braiding processes involv¬ 
ing a red base loop. We can find the braiding statistics 
associated with a blue base loop, by simply switching the 
roles of “red” and “blue”: 


°lr,b = n + n(q r + q' r ), 

e lb,b =tt (Qb+q'b), 

°br,b--^+A ( lr+q , b)- (76) 

Putting together equations dm4m<!isD and , we 
have found all the three-loop braiding statistics of the 
Hi model. (Actually, one quantity that we have not com¬ 
puted is the exchange statistics of the loops. However, 
this quantity is not necessary for our purposes, since the 
mutual statistics computed above is sufficient to distin¬ 
guish the two models). 

Before concluding this section, we would like to men¬ 
tion that the three-loop statistics for Hi appears to 
be identical to the three-loop statistics of one of the 
Z 2 x Z 2 Dijkgraaf-Witten models2&2&. To see this, we 
follow the notation of Ref. S3 and M and we summarize 
the three-loop statistics for Hi in terms of the quantity 
Qij^k = where the indices i, j, k G {r, fr}, and Sij,k 

is defined modulo 2i r. Translating equations ®,dZD, 
and into the 0 variables gives: 

066, r — 0 rr,r — 0rr,6 — 0 66,6 — 0 

0r6,r — 06r,6 — 7T, (77) 

If we compare this data to the results of Ref. Hi we 
can see that it matches the three-loop statistics for the 
Z 2 x Z 2 Dijkgraaf-Witten model labeled by (pi = l,p 2 = 
1). Based on this fact, we conjecture that Hi belongs to 
the same phase as this Dijkgraaf-Witten models 


3. Comparing the three-loop statistics in the two models 

With the above results, we now show that the Hq and 
Hi models have distinct three-loop statistics. To see an 


example of a difference between the two models, consider 
the formula for 0^ br (T75|) . This formula implies that, in 
the Hi model, if we braid a red loop around a blue loop, 
while both are linked to a red loop the resulting phase 
is ±7r/2 depending on what type of red and blue loops 
are being braided. In contrast, in the Hq model, we can 
see from (l6l|h (l62|h (l63|) that if we braid any two loops 
around one another, the resulting phase is always 0 or 7r, 
independent of the choice of base loop or what loops are 
being braided. 

For another example of a difference, consider the for¬ 
mula for 0l b r for the case where the two blue loops 
(fr, qb,q r ), (fr, q'b, q' r ) are identical , i.e. qb = q b and q r = q' r . 
In this case, Ok r = 7T. This means that if we braid two 
identical blue loops around one another while they are 
both linked to a red loop, the resulting statistical phase 
is 7r in the Hi model. On the other hand, in the Hq 
model, we can see from (I6T1) . (l62]h (l63|) that if we braid 
any two identical loops around one another the resulting 
phase is always 0, independent of which loops or base 
loops are involved. 

From the above examples, it is clear that there is no 
way to map the loop excitations of Hq onto the loop 
excitations of Hi in such a way that the correspond¬ 
ing loops have the same statistics. Importantly, we can 
rule out all possible mappings between the two sets of 
loop excitations including those that change the “col¬ 
ors” or “charges” of the loop, e.g. that map the exci¬ 
tation (fr, 1,0) in H 0 onto (r, 0,1) in Hi. (This generality 
is important because the “colors” and “charges” of the 
loops have no physical meaning in this context except as 
a scheme for labeling excitations). We conclude that the 
two models have physically distinct three-loop statistics, 
and hence must belong to different phases. 


V. CONCLUSION 

In this paper, we have presented an explicit computa¬ 
tion of the three-loop braiding statistics of two spin mod¬ 
els, Hq,Hi. The key step in our analysis was our con¬ 
struction of membrane operators that create and move 
loop excitations. With the help of these operators, we 
were able to implement the three-loop braiding process 
on the lattice and find the associated statistical phase 
in each of the models. While technically complicated, 
this membrane operator approach has the advantage that 
it is more direct than previous calculations based on 
dimensional reduction^ or modular transformations on 
a 3D toruS i 15 d 6 An additional feature of our results is 
that they provide a concrete demonstration of the util¬ 
ity of three-loop braiding statistics for distinguishing 3D 
gapped phases: indeed, we have shown that the two mod¬ 
els Hq and Hi share the same particle exchange statistics 
and the same particle-loop and loop-loop braiding statis¬ 
tics. The only way to see that these models belong to 
distinct phases is to examine their three-loop braiding 
statistics. 
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The discussion in this paper raises an important ques¬ 
tion, namely whether three-loop braiding statistics, to¬ 
gether with particle exchange statistics and particle-loop 
braiding statistics, provides enough data to uniquely dis¬ 
tinguish all 3D gapped phases. The existence of the 3D 
cubic code models suggests that the answer to this ques¬ 
tion may be ‘no’ in general. Indeed, the cubic code model 
does not support deconfined particle-like and loop-like 
excitations like the models studied in this paper, so it is 
not clear whether one can define particle or loop braid¬ 
ing statistics for this system. In light of this example, a 
more natural question may be whether the above braid¬ 
ing statistics data is sufficient to distinguish an appro¬ 
priate subset of 3D gapped phases, such as the “gapped 
quantum liquids” defined in Ref. |40|. As far as we know, 
this is an open question: we are not aware of any coun¬ 
terexamples or arguments one way or the other. 

One direction for future work would be to investigate 
the implications of our results for symmetry protected 
topological (SPT) phases. Indeed, as discussed in the in¬ 
troduction, both Hq and Hi can be thought of as Z2 x Z2 
gauge theories obtained by gauging the Z2 x Z2 symme¬ 
try of two different spin models. The Hi model comes 
from a spin model^i belonging to a nontrivial SPT phase, 
while Ho comes from a spin model in a trivial SPT phase. 
Following the same approach as in Ref. IS it should be 
possible to derive bulk and surface properties of the two 
SPT phases from the braiding statistics in the associated 
gauge theories Ho, Hi. Results of this kind will be dis¬ 
cussed in a separate publication. 35 

Another direction would be to construct exactly solu¬ 
ble lattice models that can realize more general types of 
3D gapped phases. For example, since the two models 
Ho, Hi are built out of two species of intersecting mem¬ 
branes, it may be possible to build more general models 
by considering more species of membranes or by allow¬ 
ing membranes to branch — in analogy with the string- 
net models of Ref. 0. Given our suspicion that Ho, Hi 
belong to the same phase as Z2 x Z2 Dijkgraaf-Witten 
models, a key question is whether such a construction 
can realize any phases that cannot be realized by pre¬ 
viously known exactly soluble models, such as the 3D 
Dijkgraaf-Witten model s 20 ! 36 or 3D string-net models^. 


ACKNOWLEDGEMENTS 

We thank Chenjie Wang for useful discussions. This 
work is supported by the Alfred P. Sloan foundation and 
NSF DMR-1254741. 


Appendix A: Proving the identities (| 2411 26 II 2 7|) 

In this section, we prove a few identities involving the 
operators A^A^ F>°, F?°, B\, B\ and the wave functions 
|To) and |Ti). We used these identities when we solved 
the two models, Ho and Hi. The first set of identities 




FIG. 13. An example of the 2D identity (1A8D for the case 
of loops living on the honeycomb lattice. Panel (a) shows a 
loop configuration A, while panel (b) shows the corresponding 
configuration X + p. In this example, iVi 00 p(A) = 2 and 
Ai oop (A + p) = 1. Also n p (X) = 4, since four of the six 
legs adjacent to the plaquette p are occupied by strings. The 
identity (1A8D holds for this example, since (—l) 1-2 = (—1 )-i 4 . 

state that |To) and |Ti) are eigenstates of the Ai,Aj 
operators: 

A^o) = A^o) = |*o>, 

4|*i> = = l*i>- ( A1 ) 

Similarly, the second set of identities state that |\l/o) is an 
eigenstate of BB® and j^i) is an eigenstate of B* , Bg : 

B c °|*o> = B c °|*o> = |*o>, (A2) 

B C 1 |^ 1 ) = B C 1 |^ 1 ) = |^ 1 ). (A3) 

The last identity states that the B \, B\ operators com¬ 
mute with one another: 

[BlBt,] = [BlBl] = [BlB\\ = 0- (A4) 

We begin by deriving Eqs. (EU. These relations are 
obvious since |To) and |Ti) are linear superpositions of 
closed membrane configurations |X&, A r ), which by def¬ 
inition obey Ai\Xb,X r ) = Aj\Xb,X r ) = \Xb,X r ). The 
relation (1A2|) is also easy to prove: first we note that 

B° c \* 0 ) = Y B° c \X b ,X r ) 

closed Xb,X r 

= Y I Xb + C,X r ) 

closed Xb,X r 

= Y I X iX r ) 

closed X' h ,X r 

= |^o) (A5) 

where we made the change variables X' h = X^ + c in the 
third line. Then using the fact that F>° = ^(1 + S°), we 
immediately derive the required relation F?°|To) = |Tq). 
The same argument works for F>°. 

The relation (|A3|) requires a little more work. First, 
we note that 

Y (-l) Ng(Xb ' Xr) Bl\X b ,X r ) (A6) 

closed Xb,X r 

= Y (-l) NgiXb ’ Xr \-l) m H n c\X b + c,X r ) 

closed Xb,X r 
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where N g is the number of green loops in the configu¬ 
ration X^X r and m c ,n c are integer valued functions of 
Xb, X r which are defined in section HID 21 To proceed 
further, we use the following identity: 

(_1 \N g (X b +c,X r )-N g (X h ,X r ) = ^_pm c (X 6 ,X r )^ c (X b ,Xr) i 

(A7) 

Here, we explicitly show the dependence of m c and n c on 
Xb,X r , for clarity. Equation (IX7l) says that when we add 
a blue cube, i.e. X 5 —» X^ + c, the resulting change in the 
parity of the number of green loops N g can be computed 
in terms of m c and n c . This relation is nontrivial because 
(_l)N g depends on the global properties of the membrane 
configuration X&,X r , while m c and n c only depend on 
local properties of X\>, X r in the neighborhood of the cube 


To derive the above identity (1X71) . we first recall an 
analogous, but simpler identity involving 2D loop models. 
The 2D identity applies to 2D loop configurations X that 
live on the links of a 2 D planar trivalent lattice. It states 
that if we add a loop around the boundary of a single 
plaquette p to form a new loop configuration X +p, then 

(_ 1 )M oop N\ 00 p(x) _ ^ ^) . (A. 8 ) 


where N\ oop (X ) and N\ oop (X+p ) are the number of loops 
in the configurations X and X +p, and where n p (X) is 
the number of occupied legs of the plaquette p. (See 
Fig. [13] for an example). The 2 D identity (1A8[) is closely 
related to the 2D exactly soluble double semion model of 
Ref. 0, and can be established using the local rules that 
define that model. 

To derive the 3D identity (|A7 [) from its 2 D cousin (|A 8 |L 
it is helpful to think about the configuration of green 
loops in A^,X r , for a fixed red membrane configuration 
X r . These green loops he on the surfaces of the red mem¬ 
branes X r . When we add a blue cube, this effectively 
adds a collection of green loops along all the red mem¬ 
branes that intersect the cube. Adding these green loops 
is similar to adding the loop p in the 2D identity (IA 8 [L so 
when we add these loops, the change in the parity of the 
number of green loops can be calculated by multiplying 
the factors on the right hand side of (|A 8 [) over all the ad¬ 
ditional green loops. Taking the product of these factors 
gives the 3D identity (1A7[) . 

Now that we have established the identity (IA7l) , we 
can complete our derivation of (1A3[) : substituting (1A7[) 
into (1A6[) , we obtain 


Scl*i> 


{-l) N ^ x » +c ' x ^\X b + c,X r ) 

closed X b ,X r 

closed X^,X r 

l*i) (A9) 


where we made the change of variables X' h = X^ + c 
in the second line. Finally, using the fact that E>1 = 
|(1 + B\)P C and P c |Ti) = |4/i), we derive the required 


relation B\ |Ti) = |Ti). The same argument works for 

Bb 

All that remains is the third relation (|A4|) . In principle, 
one could establish (1A4[) by writing out the explicit form 
of B\,B\ and calculating their commutators. However 
such a calculation would be tedious and not particularly 
illuminating. Therefore we use a different approach to 
show (1A4[) . We make use of two properties of the opera¬ 
tors B\. The first property is the relation 

2£|*i> = |*i) (A10) 

which we just established above. The second property is 
that, for any closed membrane configuration (X^,X r ), 

Bl\X b ,X r ) = e id ^ x ^\X b + c,X r ) (All) 

where e l()(Xh ' Xr ’ c> is a phase factor that depends on 
W,,X r ,c. The latter property follows immediately from 
the definition of B 

Using the above two properties (1A10[) and (1A11[) , we 
will now show that the two operators B\,B\, commute 
with one another when acting on closed membrane states. 
The argument is simple. From property (1A11[) , we com¬ 
pute 

\X b , X r ) = B\, ■ e iS ( Xb,Xr ’ c }\X b + c, X r ) (A12) 

= e iO{X»+c,Xr,c') . e ie(X b ,X r ,c) | X b + C + d , X r ) 

On the other hand, if we reverse the order we find 
B\B\,\X b ,X r ) = 

e i0(X b + c',X r ,c) . e i0(X b ,X r ,c') \x b + c + c ' 5 x r )(A13) 

We need to show that the expressions on the right hand 
sides of (1A12|) and (|A13|) are equal. To this end, we note 
that property (1A10|) implies that 


(X b ,X r \B 1 c \V 1 ) = (X b ,X r \* 1 ). (A14) 

Then, using property (1A11I) . we can write this as 

e -«>(X b ,X r ,c)( Xb + CfXr 1 ^) = (X^Xrl*!) (A15) 

so that 


i6(X b ,X r ,c ) _ (Xb + C,X r \^i) 

(X^Xrl*!) 


(A16) 


We now substitute (|A16|) into equation (1A12[) . The result 
is 


B\,B\\X b , X r ) = I X b +C + c', X r ) • 


(X 6 + c + c',X P |* 1 ) 

(X b ,X r |^!) 


Likewise, if we substitute (IA16I) into equation (1A13I) . we 
obtain the same expression: 


B\B\,\X b ,X r ) 


\X b + c + d, X r ) 


{Xh+C + d.Xr^dl 

(X^Xr^x) ’ 





























21 


implying that Bl and B\, commute when acting on closed 
membrane states. 

Now that we know that B\ and B\, commute when 
acting on closed membrane states, we can immediately 
deduce that the combinations B\P C and B\,P C > commute 
when acting on arbitrary membrane states, since the op¬ 
erators P c and P c f project onto states that are closed in 
the neighborhood of c and d. It then follows that the 
operators B\ and B\, commute with one another, since 
B\ is a linear combination of B\P C and P c . 

We have shown that = 0 for any cubes c, d 

in the cubic lattice. This establishes the first identity in 
Eq. (jAU). The other two identities in Eq. dm can be 
proven in exactly the same way. 


Appendix B: Showing that the cylinder operators 
create braiding eigenstates 


number (in fact, a phase factor). We now translate this 
claim into a mathematical equation. 

Consider, for example, braiding another blue loop 
around one of these loops. This braiding process can be 
implemented by applying a blue torus operator M^(S') 
where the torus S' encircles the loop that we wish to 
braid around (see Fig. [12] for a similar geometry). The 
fact that the blue loop is an eigenstate of this braiding 
process is equivalent to the equation 

■ (Mt(S)\V ex )) = (const.) • M£(S)\V ex ). (B2) 

Our task is to show that if f b obeys m then the 
operators obey equation (IB21) . To establish this re¬ 
sult, we first rewrite (1B21 in a more convenient form. 
We begin by recalling that since S' is a torus opera¬ 
tor, it does not create any excitations — in other words, 
Ml {S') |T ex ) oc | T ea ,). Using this fact, the above equa¬ 
tion can be written as 


In this section we show that the blue cylinder operator 
Ml(S) creates loop excitations which are eigenstates of 
braiding measurements. Our argument consists of two 
parts. In the first part, we show that as long as f b obeys 
condition ([46]) . then the loop excitations are guaranteed 
to be eigenstates of braiding measurements. In the sec¬ 
ond part, we check that f b does in fact obey ([46]) . The 
second part can also be thought of as a derivation of 
Table [J and Table mi since we will see that the values in 
those tables are completely fixed by the requirement that 
f b obey equation (146]) . 


1. Connection between equation (|46|) and braiding 
eigenstates 


Mt(S')Mt(S)\V ex ) = (const.) • M^(S)M^S')\V ex ). 

Next, multiplying both sides by (X b ,X r \ gives 

{X b ,X r \Ml(S')Ml(S)\y ex ) = 

(const.) • {X b ,X r \Mt{S)Mi(S')\y ex ). 

Using the definition of M^ 1 , we rewrite the matrix ele¬ 
ments on the left and right hand sides as: 

(X b ,X r \Mt{S')Mt{S)\* ex ) = ft(X b ,X r ,S') 

• ft(X b + 5', X r , S ) • * ex (X b +S + S', X r ), 
(X b ,X r \Mt{S)Mt(S')\9 ex ) = f;(X b ,X r ,S) 

• ft{X b + 5, X r , S') • y ex (X b + S + S',X r ). 


In this section we show that if f b obeys condition (146]) 
then the blue cylinder operator is guaranteed to cre¬ 
ate loop excitations which are braiding eigenstates. We 
reprint equation (146]) below for convenience: 


Therefore, what we need to prove is 

fb(X b , X r , 5') • f b {X b + 5', X r , S) = 

(const.) . f b {X b , X r , 5) • f h {X b + S, X r , S') 


f b {X b ,X r ,Si U S 2 ) = f h {X h ,X r ,S i) . f h {X h ,X r ,S 2 ). 

(Bl) 

Here Si and S 2 can be any two cylinders that share 
a common boundary while {X b ,X r ) can be any mem¬ 
brane configuration whose intersection with the common 
boundary obeys the relevant cylinder operator boundary 
condition. 

To begin, it is helpful to consider an example. Let 
| T ea ,) be an excited state with a single red loop. We will 
call this red loop the “base loop.” Let M^(5) be a blue 
cylinder operator that is linked with the base loop. Con¬ 
sider the state Ml(S) |T ex ). This state contains two blue 
loops at the two ends of the cylinder S', both of which 
are linked with the red base loop. Now, what we want 
to show is that these blue loops are eigenstates of braid¬ 
ing measurements. That is, we want to show that if we 
braid any other loop around one of these loops, the sys¬ 
tem returns to its original state multiplied by a complex 


or equivalently 


f b (X b + S',X r ,S) f b {X b ,Xr,S') 
f b {X b ,Xr,S) ’ f h (X h + S, Xr, S') 


(const.). (B3) 


We now prove (1B3D . Focusing on the first ratio in 
(1B3D , we note that the numerator and denominator can 
be represented graphically as 


f b {X h + S',X r ,S) = f b 
f b {X b ,Xr,S)=f b 



Here the thin vertical blue region corresponds to the in¬ 
tersection S' fl 5, while the gray areas on the left and 
on the right are meant to denote some arbitrary pictures 
corresponding to X b D S and X r H S. (We can assume 
without loss of generality that the pictures are of the 
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above simple form, since we can use the path indepen¬ 
dence property, discussed below Eqs. ([36]) . to deform S 
and S' as we wish). Next we use equation m twice 
to break up the cylinder S into three smaller cylinders, 
giving us 



Taking the ratio of these expressions gives 


ones below: 



Next, we use the constraint equations ([36]) to reduce the 
pictures in the second and third lines to basic pictures: 


f b (X b + S',X r ,S) 
f b (X b ,X r ,S) 


fb 

h 





\ 




/ 



\ 



/ 


In particular, we see that independent 

of X\),X r . In exactly the same way, one can show that 
ff(x^+sx S S ') independent of Xb,X r . This proves 
equation (|B3D : both terms on the left hand side are inde¬ 
pendent of Xfc, X r , so their product must be independent 
of X\>, X r as well. 


So far, we have shown that m implies the braiding 
eigenstate property in one case. The case we considered 
was a blue loop linked to a red base loop. We showed 
that the blue loop is a braiding eigenstate with respect 
to braiding another blue loop around it. To prove the 
general claim, we need to establish the same braiding 
eigenstate property for all other cases, that is all other 
combinations of red and blue loops. We will not discuss 
the other cases here, but it should be clear that the above 
argument is not specific to the case considered above and 
applies equally well to the other cases. 


2. Checking that fb obeys equation (|46|) 


In this section, we check that fo obeys equation mi) 
— or equivalently m- To begin, we note that it suffices 
to check m for the case where the intersections X & D Si 
and X r n Si look like one of the four basic pictures shown 
in Table [T] and Table an The reason that we only need 
to consider this case is that the constraint equations HMD 
that define fb are local, so if m holds for basic pictures, 
then it automatically holds for general Xb,X r . 

Let’s start with the simplest case: blue cylinder oper¬ 
ators that create unlinked loops. In this case, the four 
basic pictures are those shown in Table [U Specializing to 
these pictures, equation m reduces to 4 2 = 16 different 
relations that we need to check. Many of these equations 
are redundant, so we only write the four independent 


fb 

fb 



Putting this together, we have 
fb (\ \) * fb 


* (□□)•* 

fb 


fb 


• fb 

• fb 


)=//' 
)- a 
)=/. 
)=// 


Letting 


h 




yen) 

] 


fj n 

vJ 

h \ JJ 


= D, 


we can rewrite these equations as 

A 2 = A, 
B 2 = A, 
C 2 = A, 
B ■ C = D. 


(B4) 


Now let us compare with the values for A, B,C,D given 
in Table H 

A = 1, B = -e inqr , C = e inq \ D =-e^ qr+qh) 

(B5) 

where q r = 0 or 1. First, we can see that the above 
expressions (|B5j) do in fact obey equations (|B4j) . Thus, 
we have successfully verified 4HU> for the case of an un¬ 
linked blue cylinder operator. In fact, we can see that 
(|B5|) is actually the most general possible solution to the 
equations (lB4l) . Thus, the above calculation can also be 
regarded as a derivation of Table [H 
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Next, we consider the case of a blue cylinder operator 
linked with a red base loop. We proceed in exactly the 
same way as in the unlinked case. First, we note that 
equation (ED reduces to 4 2 = 16 different relations, of 
which 4 are independent: 


Letting 


fb 

fb 


= A 
= c, 



= B, 
= D, 



Next, we use the constraint equations dnni) to derive: 


fb 

f b 


- ~ h (B) • 


Putting this together gives 


fb 

fb 

fb 

fb 



• fb 
•fb 


* fb 

• fb 


VQ 





= - fb 
= ~fb 



we arrive at the following algebraic equations: 

A 2 = A, 

B 2 = -A, 

C 2 = -A, 

B ■ C = D. (B6) 


Let us compare with the values of A, B, C, D given in the 
top panel of Table [Tl] 

A = 1, B = i ■ e i7rQr , C = i ■ e i7Tqb , D = -e^+q*). 

(B7) 

Again, we can see that the above expressions ED obey 
(|B6[) . Thus, we have proven equation ED for the case of 
a blue cylinder operator linked to a red base loop. We can 
also see that the above expressions are the most general 
solutions to these equations. Thus, our calculation can 
also be thought of as a derivation of the top panel of 

Table HI 

Finally, we consider the case of a blue cylinder operator 
linked with a blue base loop. Proceeding in the same 
away, we set 


fb 

fb 


Vet 


= A 

= c. 



and we derive the algebraic equations 


A 2 = A, 
B 2 = A, 

C 2 = A 
B-C = D. 


(B8) 


Again, it is easy to check that the values of A, B,C,D 
given in the bottom panel of Table U obey equations 
(|B8 [) . This establishes equation (jBlj ) for the case of a 
blue cylinder operator linked with a blue base loop. 
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